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ABSTRACT 



The Supersymmetric SO(10) GUT based on the 210 © 10 © 120 © 126 © 126 
Higgs system and only spontaneous CP violation has 24 real superpotential param- 
eters. It is compatible with all current low energy data and yet has enough free- 
dom to accommodate the still unknown Leptonic CP violation and neutrino mass 
scale parameters. It is thus the New Minimal Supersymmetric Grand Unified The- 
ory (NMSGUT). We compute the superheavy mass spectrum, couplings and associ- 
ated gauge evolution threshold effects and determine the regions of the parameter 
space compatible with perturbative unification, baryon stability and Weinberg angle 
while permitting spontaneous CKM CP violation. We show that the structural re- 
quirements of spontaneous CP violation always tend to raise the unification scale Mx 
and with it the mass of all superheavy multiplets and thus to suppress Baryon de- 
cay. These scales can be raised even close to the Planck scale so that the NMSGUT 
gauge Landau pole becomes innocuous. This scenario lends support to our earlier 
speculations regarding GUT dynamical symmetry breaking and induced gravity via 
UV condensation at this Landau pole. Fermion mass formulae and the effective su- 
perpotential for Baryon violating d = 5 processes - which contains novel channels due 
to the 120 Higgs - are derived to prepare for the next round of detailed analysis of 
the phenomenology of this model. 



1 



1 Introduction 



The discovery of neutrino mass was both preceded bypQ and itself provokedjH [31 HI 
13 El [7] intensive investigation of unifying theories that naturally incorporate super- 
symmetry and the seesaw mechanisms [8 J : in particular models with the Left-Right 
gauge group as a part of the gauge symmetry and B — L broken at a high scale 
and R/M-parity preserved to low energies [Tj. The close contiguity of the seesaw 
scale and the Grand Unified scale pushed SO(10) GUTs, which are the natural GUT 
home of both Type I and Type II seesaw mechanisms, but were long relegated as 
baroque cousins of the -seemingly- more elegant minimal SU(5) GUT, into centre 
stage. The understanding that the Susy SO(10) based on the 210 © 10 © 126 © 126 
Higgs system proposed [TOl [TT| long ago was the best candidate for the Minimal Super- 
symmetric GUT (MSGUT) crystallized after the demonstration of its minimality on 
parameter counting grounds and an elegant reduction of its spontaneous symmetry 
breaking problem to a single cubic equation with just one unknown parameter |12j. 
Careful computations of the symmetry breaking and mass spectra [101 [HI [12] became 
available [31 HI El E] MSGUT. These theories naturally maintain a structural distinc- 
tion between Higgs and matter fields and therefore naturally preserve R-parity down 
to low energies [131 [TJ [2]. 

The initial euphoria[7j that the version utilizing only 10, 126 Higgs representations 
might prove sufficient [9 J to fit all low energy fermion data in an elegant and predictive 
way ran aground when the successful generic fits of fermion mass data were shown to 
be unrealizable [Til [T5] in the context of the actual Seesaw mechanisms available in 
the MSGUT. Thereafter the GUT based on the 210 © 10 © 120 © 126 © 126 Higgs 
system has emerged |14l [151 [TBI [T71 [18] as the New Minimal Supersymmetric GUT 
(NMSGUT) capable of fitting all the known fermion mass and mixing data. Re- 
markably this compatibility has been demonstrated even when all the superpotential 
parameters are real and with 3 matter yukawas set to zero and the 126 vev fixed in 
an ad-hoc way |18] . Thus at least 4 free parameters are still available to fit the still 
unknown values of the leptonic CP violation phases, neutrino mass scale and 0f 3 when 
they become available. The theory with the 120 thus escapes the contradiction with 
the composition of the light Higgs doublets which was demonstrated [HI [15] to occur 
in the MSGUT without it. Thus this actual or New MSGUT calls for and deserves 
the same detailed analysis of its superheavy threshold effects, fermion mass fit com- 
patibility and exotic effect effective superpotential that we earlier provided for the 
theory without the 120 [IDl 151 151 H3| H5j which we previously called the MSGUTp] 
but whose claim on that name is now tenuous and faded. 

Even after it fits the fermion data the NMSGUT must still face the challenge 
posed by the non observation of Baryon decay in spite of the lore that minimal 
Susy GUTs imply baryon decay rates via d — 5 operators which are higher than 
the current experimental upper bounds[2Tj. Moreover all the analysis of MSGUTs 
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has so far been done in the framework of softly broken global supersymmetry with 
no attempt to incorporate it's dynamical or spontaneous breaking. Furthermore the 
NMSGUT suffers [22] from a Landau pole at Ax lying just above the perturbative 
unification scale Mx due to the huge SO(10) gauge beta functions implied by the 
large representations used. 

Last week, as this paper was being finalized, a very interesting paper [23] appeared 
which gave a lucid discussion of these essential and perennial issues [21] and suggested 
that they could be addressed simultaneously by considering metastable Susy breaking 
vacua in the context of N=l Supergravity with non-renormalizable operators included 
in the matter superpotential to allow escape of the unification scale and d = 5 baryon 
decay amplitudes from the tight confines [2Tj of the minimal SU(5) SGUT. The au- 
thors argued that the MSSM singlets of the MSGUT itself could play the role of 
the supersymmetry breaking field(s) of the Polonyi or other hidden sector and that 
this melded optimally with a Mx ~ Mpi anc k solution to the d = 5 baryon decay 
rate difficulty of Susy GUTs. For reasons we shall explain below we do not find a 
non-renormalizable superpotential cogent or attractive inasmuch as it puts paid to 
the predictivity engendered by renormalizability. It seems moreover that it has been 
invoked more to provide additional parameters that may be fine tuned to satisfy the 
needs of the attractive phenomenological scenario suggested [23j. It is natural to en- 
quire whether the same program might be naturally fulfilled in a realistic theory with 
only a renormalizable superpotential such as the MSGUT for which we had computed 
the threshold corrections alreadyjHJ [HI [15] and the NMSGUT which we were in the 
process of reporting. In fact such a raising of the unification scale was observed by us 
in considerable regions of the parameter space already [5], [HI [15]. Since however the 
focus had been on raising neutrino masses by lowering the right handed neutrino mass 
we had not yet questioned whether the new found [HI [17] leeway in fitting neutrino 
masses via Type I seesaw with highly suppressed 126 yukawa couplings [151 US ELY] 
would allow us to exploit the domination of the raised unification scale regions in the 
viable parameter space of the MSGUT and NMSGUT. We therefore remedied this 
defect and report on the very encouraging indications in detail below. We note that a 
novel type of "intermediate" scale Dirac supermultiplet supermultiplet transforming 
as [8, 2, ±1] with respect to the SM gauge group - which couples to quarks via the 
126 yukawa - becomes lighter as the unification scale is raised. It is the effects of 
this multiplet that will alter the RG equations sufficiently to raise Ac consistently. 
Unlike similar multiplets used in[23] for the same purpose, however, the mass of this 
multiplet does not seem to fall below 10 15 GeV and that too only when Mx and 
have long become unphysical. 

We have speculated [25J that the vice of Asymptotic Strength(AS) of (N)MSGUTs 
may be turned to good account to construct a theory of dynamical calculable GUT 
symmetry breaking using an extension of the techniques for dealing with strongly cou- 
pled supersymmetric theories[26], and even that AS Susy GUTs(ASSGUTs), which 
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determine by their RG flow their own physical UV cutoffs, escape the objections to 
the induced gravity program that halted its development in the 80's[27]. This induced 
gravity program would gain plausibility if the Planck scale and Ax coincided. Such 
a coincidence is an obvious consequence of raising the perturbative unification scale 
to just below the Planck scale. We shall see that the requirement of viable unifi- 
cation in the presence of high scale threshold corrections in the NMSGUT leads us 
almost inevitably to regions of the parameter space where Mx, together with baryon 
decay mediating triplet masses, are raised above 10 16 GeV and A x approaches closer 

tO M Pl anc k ■ 

The computation of mass spectra, couplings and gauge threshold effects in the 
NMSGUT and a partial analysis of their implications is the subject of this paper. 
A calculation of the chiral mass spectra using different conventions and method is 
already available [31]. However, since our methods also yield all couplings which we 
need for an evaluation of the phenomenological implications, we present all results 
consistently using our own conventions. 

In Section 2 we recapitulate our notation and the basics of such models. Detailed 
accounts of our techniques have already been given earlier [3], [5], [H]. In Section 3 
and Appendix A we give the mass spectra and in Appendix B we describe an SU(5) 
reassembly crosscheck of the spectra we obtain. In Section 4 we discuss how the 
threshold effects calculated using these spectra determine regions where perturba- 
tive unification is viable and sub-regions of these where Baryon decay is naturally 
suppressed. We shall see how both in the general case of complex superpotential pa- 
rameters and when the CKM CP violation is spontaneous one is generically led, in the 
latter case inescapably, towards a raised unification scale: thus leading to a potential 
resolution of the basic problems of Susy GUTs discussed above, without any con- 
trived cancellations. In Section 5 we give the fermion mass formulae in the presence 
of the 120 plet using analytic expressions for the null eigenvectors (after fine tuning 
to keep one pair of doublets light) of the 6x6 Higgs doublet ([1, 2, ±1]) mass matrix 
(Appendix C). In Section 6 we integrate out the heavy triplets that mediate Baryon 
decay via d = 5 operators and give the resultant effective superpotential in terms of 
the matter superfields of the effective MSSM. We conclude with a brief discussion of 
issues and on-going/future work in Section 7. 

2 The New Minimal Susy GUT 
2.1 MSGUT couplings, vevs and masses 

The original MSGUT [TOl [TT| [9[ [12] was the renormalizable globally supersymmetric 
50(10) GUT whose Higgs chiral supermultiplets consist of AM type totally antisym- 
metric tensors : 2 10 126(£^Hm), 126(£jjM m )(i, j = 1...10) which break the 
GUT symmetry to the MSSM, together with Fermion mass (FM) Higgs lO-plet(Hj). 
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The 126 plays a dual or AM-FM role since it also enables the generation of realistic 
charged fermion and neutrino masses and mixings (via the Type I and/or Type II 
Seesaw mechanisms); three 16-plets ^a(^4 = 1, 2, 3) contain the matter including the 
three conjugate neutrinos ipf;)- 

The superpotential (see[T2l 131 HI 15] for comprehensive details ) contains the mass 
parameters 

m : 210 2 ; M : 126 • 126; M H : 10 2 (1) 
and trilinear couplings 

A : 210 3 ; 7] : 210 • 126 126; 7 © 7 : 10 • 210 • (126 © 126) (2) 

In addition one has two complex symmetric matrices Hab, Jab of Yukawa couplings 
of the the 10, 126 Higgs multiplets to the 16.16 matter bilinears. The U(3) ambiguity 
due to 5*0(10) 'flavour' redefinitions can be used to remove 9 of the 24 real parameters 
in /, h. In addition rephasing of the remaining 4 fields H, S, S removes 4 phases 
from the 14 parameters in m, M, Mjj, A, rj, 7, 7 leaving 25 superpotential parameters 
to begin with. Strictly speaking, since a fine tuning to keep one pair of doublets light 
is an intrinsic part of the MSGUT scenario, an additional parameter (say M# ) may 
be considered as fixed so that there are actually 24 free superpotential parameters. 
In addition the electroweak scale vev vw, and tan/3 are relevant external parameters 
for the light fermion spectrum determined by the GUT yukawa structures. 

The GUT scale vevs that break the gauge symmetry down to the SM symmetry 
(in the notation of [3]) are[10l [H] 

(3) 

(4) 
(5) 
(6) 
(7) 

The vanishing of the D-terms of the SO(10) gauge sector potential imposes only the 
condition \a\ = |cr|. Except for the simpler cases corresponding to enhanced unbroken 
gauge symmetry (577(5) x U(l), 577(5), G 3>2) 2,b-l, Gz,2,r,b-l etc) [12111], this system 
of equations is essentially cubic and can be reduced to a single equation [T2] for a 
variable x = —Xu/m, in terms of which the vevs a, u,p, a, a are specified : 

C(x, = 8s 3 - 15x 2 + Ux - 3 + f (1 - xf = (8) 

where £ = — . Then the dimensionless vevs in units of (m/A) are Q = —x [12] and 

(x 2 + 2x-l) x(5x 2 -l) ~ 2 Xx(l - 3x)(l + a; 2 ) 

P = — Za2~ 5 acy =- 7T— 7a ( 9 ) 
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This exhibits the crucial importance of the parameters £, x. Note that one can 
trade [121 111 US] the parameter £ for x with advantage using the equation (jSJ) since £ 
is uniquely fixed by x. By a survey of the behaviour of the theory as a function of 
the complex parameter x we thus cover the behaviour of the three different solutions 
possible for each complex value of £. However the £ parametrization is crucial to 
study the spontaneous CP violation case that is a principal focus of this paper. 

2.1.1 Digression on Supergravity Potential 

So far the vacuum of MSGUTs has been determined by extremizing the global susy 
potential by satisfying the F and D conditions [TQl [HI 12] ■ This was justified since 
the supergravity corrections to F would be ~ Mjr/Mp « M\, if M x ~ \Q lG GeV 
as is customary in Susy GUTs. However in view of the proposal |23] to raise the 
GUT scales to the vicinity of the Planck mass, by virtue of threshold corrections 
that modify the MSSM unification value, Supergravity corrections to the potential 
would surely be important. If moreover as proposed in |23j one could employ the 
MSSM singlets in the GUT fields themselves to break supersymmetry and to tune 
the superheavy contributions to the cosmological constant to zero, the minimization 
of the full supergravity potential becomes essential. 

In the MSGUT and NMSGUT the there are 5 MSSM singlet chiral fields namely 
p(l, 1, 1)210, a(15, 1, 1)210, ^(15, 1, 3)210, er(10, 1, 3) i2 6, cr(10, 1, 3) m (in terms of the 
SU(4) x SU(2) L x SU(2) R labels of the submultiplets they lie in). The MSSM 
singlet superpotential takes the form 

W(a,p,u,a,a) = m(p 2 + 3a 2 + 600 2 ) + 2A(a 3 + 3pu 2 + 6au 2 ) 

+ (M + ri(p + 3a-6u))aa + W + W dyn (10) 

where we have included a constant term Wq and a dynamical superpotential [25] 
function of chiral SO (10) invariants formed from the the singlet fields to account for 
gaugino condensation and non trivial chiral moduli when the theory condenses in 
the UV. Such a superpotential was used to demonstrate the feasibility of (calculable) 
dynamical symmetry breaking in an asymptotically strong 377(2) toy model[25]. This 
superpotential is highly constrained by the holomorphy of the Supersymmetric gauge 
theory and can in principle be determined by arguments similar to those we used for 
the toy model [25]. We propose, in good time, to do the same for the MSGUT but 
for the moment we shall ignore it. No other renormalization of the superpotential is 
expected due to holomorhy[26j. 

We have deliberately not inserted any other SO(10) invariant non-renormalizable 
gravity induced terms in the superpotential since we expect that gravity itself may 
be induced by the UV condensation of the GUT. In the aftermath of the AdS-CFT 
conjecture the expectation of an effective induced gravity as a consequence of gauge 
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condensation has become almost passe [28J . However here no exotic extra dimensional 
context is involved but only the one context where supersymmetry can confidently be 
expected to be obeyed and a GUT that has passed a series of consistency tests and 
challenges to its accommodative abilities without falling into a disordered state of ad- 
hoc invocation of special mechanisms invoked to solve one-off problems. Predictability 
may well be retained since the induced gravity should have a potential of the N = 1 
supergravity form which is determined by the same unrenormalized (but dynamically 
corrected for gauge condensation) superpotential and possibly additional corrections 
to the kinetic terms by modification of the Kahler potential only [25]. If one was 
forced to include all possible or even some part of the possible gauge invariant but 
non-renormalizable superpotential operators we would lose all hope of a predictive or 
quantitative understanding that was not merely a multi-parameter re-formulation of 
our ignorance and hope. The scenario we advocate provides at least the skeleton of 
a rationale for the absence of all such arbitrary corrections due to the known power 
of holomorphy in protecting the superpotential. 

The supergravity potential takes the form [30] (k = (Mp)" 1 ; M* P = M P /\f&ix = 
2 x 10 l8 GeV) 



|2 



ft = -^(rtft + k^" 1 )!^ 1 (ii) 

K 4 



Where f a b is the gauge kinetic function and the notation is standard|33j. Clearly 
the analysis of the full problem is quite daunting. Let us begin by writing down the 
potential in the simplified case where we take the kinetic terms to be canonical and ig- 
nore the dynamical superpotential and gaugino bilinears in T . Then the supergravity 
F terms take the form [K = \p\ 2 + 3\a\ 2 + 6\u\ 2 + l£litl£li) 



T v = (2mp + QXoo 2 + r]aa + K 2 p*W)e~ K K 

T a = V3(2ma + 2A(a 2 + 2uj 2 ) + i]aa + K 2 a*W)e- K ' 2K 

= V6(2mcu + 2cu\(p + 2a) - i)oo + K 2 uj*W)e~ n2K (12) 

T a = 4=(( M + V(P + 3a - 6u)a + K 2 o*W)e~ K * K 
v2 

T s = -r{{M + r]{p + 3a-6u)a + K 2 a*W)e- K2K 
v2 

V = e « 2K (£\F i \ 2 -3K 2 \W\ 2 ) + ^J2( Da ) 2 ( 13 ) 

i Z a 
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Now the standard procedure [30] 123] can be followed : extremize dV/ d(a, p, u, a, a) = 
0, impose V\ — at the extremae and check that the mass spectrum is non tachyonic 
by appropriate conditions on the superpotential parameters. Since W is only cubic in 
the fields, the potential that must be minimized is effectively of degree 8. Although 
complicated this is definitely an improvement over a quintic [23] superpotential (with 
many unwanted relatives waiting in the the undergrowth for the door to be opened!). 
The multiplicity of singlets effectively raises the degree without introducing new pa- 
rameters. 

Assuming the D terms continue to vanish separately (only \a\ = \a\ is needed) 
and the vacuum energy at the minimum is tuned to zero, the equations which must 
be solved in the case of a cubic superpotential are (di = d/dcf) 1 ) 

didjWdjW* + tJdidjW&W* + K 2 diW(W {3) + k 2 KW)* + K 4 (j)*\W\ 2 = (14) 

Here is the part of the superpotential cubic in the fields. While a detailed 
analysis must be deferred [35] it is clear that this potential has 5 free parameters i.e 
Wq, A, rj, m, M and these could well be sufficient to obtain a solution of the desired 
type i.e with vevs < n" 1 and rri3 « Mx since the raising of gauge scale and the 
baryon decay mediating triplet will be found(see below) to be practically automatic. 
We now continue with the analysis of the theory with the global susy potential with a 
particular view to show that threshold corrections tend to raise the unification scale. 

2.1.2 Characteristics of the SSB solutions 

A knowledge of the 3 solutions Xi{^){i = 1,2,3 \C(xi(£,),£) = 0) of the cubic equa- 
tion® is necessary when considering issues such as the viability of spontaneous (GUT 
scale) CP violation in this theory. 

In Fig. 1,2 we exhibit plots of £«(£) for real £. For £ < —27.917 all three solutions 
are real. Moreover it is clear that for |se| >> l,x ~ — £/8 specifies one real branch for 
real £ as visible in Fig. 1. On the other two (complex conjugate) branches the real 
and imaginary parts of x(£) are bounded above and below( real part G (.8, 1.0) and 
imaginary part magnitude e [0, 1.1). 

Using the above vevs and the methods of [3J we calculated the complete gauge 
and chiral multiplet GUT scale spectra and couplings for the 52 different MSSM 
multiplet sets falling into 26 different MSSM multiplet types (prompting a natural 
alhpabetization of their naming convention [3]) of which 18 are unmixed while the 
other 8 types occur in multiple copies which mix. The (full details of these) spectra 
may be found in[3l [5] and equivalent results(with slightly differing conventions ) 
are presented inpEj. A related calculation with very different conventions has been 
reported in |6j. The initially controversial relation between the overlapping parts of 
these papers was discussed and resolved in [20] . 
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Re[x] 




Figure 1: Solutions of eqn.® which governs GUT ssb: Plot of Re[xi(£)] vs £ for 
i = 1,2, 3. The vertical straight line segments are "reconnection artifacts" induced 
by a switch over between real and complex solutions and vice versa. 



Among the mass matrices is the all important 4x4 Higgs doublet mass matrix 
PH E] 7i which can be diagonalized by a bi-unitary transformation [T2l HI [5]: from the 
4 pairs of Higgs doublets arising from the SO(10) fields to a new set H® 

of fields in terms of which the doublet mass terms are diagonal. 

U T HU = Diag{m { H ] ,m { H\....) 

h® = UijH® ■ h® = UijH^ (15) 

To keep one pair of these doublets light one tunes M H so that DetTC = 0. In the 
effective theory at low energies the GUT Higgs doublets h®,h® are present in the 
massless doublets in a proportion determined by the first columns of the 

matrices U, U : 

E«M X : h {i) -> ai H {1) ; a i = U <l 

h® -> ; a i = U il (16) 

The all important normalized 4-tuples a, a can be easily determined |12l 141 [T4" t [T9"j [T5] 
by solving the zero mode conditions: Tia = ; a T 7i = 0. 
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Figure 2: Solutions of eqn.(jH]) which governs GUT ssb: Plot of Im[xi{^)] vs £ for 
i — 1, 2, 3. The vertical straight line segments are again 'reconnection artifacts'. 



2.2 Additional terms introduced by the 120 

The introduction of the 120-plet Higgs representation leads to new couplings in the 
superpotential . The additional terms are 

TT1/ k p 

WnmSG — -^-^yjOijkOijk + 7yOijkH m Q 

+ 2^yjO i j] t ^] t l mn (C'El mni j + CX'lmnij) + 'AB^^Cf' 'j^ B \ x i 2 i 3 

The Yukawa coupling g^B is a complex antisymmetric 3x3 matrix. The SU{4) x 
SU{2)l x SU{2)r ( Pati-Salam ) decomposition of the 120plet is as follows : 

Otf fc (120) = O«(10, 1,1) +0^(10, 1,1) +0^(15, 2, 2) 

+ 0^(6,l,3) + 0^(6,3,l) + O ad (l,2,2) (17) 

The arbitrary phase of the 120 reduces the effective number of the extra couplings 
(m ,p, (,C,k (5 C — lr = 9 r ) and <7ab(3 c = 6 r )) so they amount to 15 additional 
parameters. Thus the relative advantage [T21 [29] with respect to SU(5) theories using 



10 



additional fields or higher dimensional operators to correct the fermion mass relations 
of the simplest SU(5) model seems weakened. However this is too hasty a conclusion. 

In fact the old MSGUT fails to fit the fermion mass data due to difficulties with 
the overall neutrino mass scale [TU [15]. An alternative scenario within the NMSGUT 
which successfully evades the problems of the old MSGUT was proposed and elabo- 
rated in [151 EE E] • m this scenario the Yukawa couplings (/ab) of the 126 are much 
smaller than those of the 10, 120. This boosts the value of the Type I seesaw masses 
(which were in any case dominant over the Type II seesaw masses but still too small) 
so that they are comfortably generically viable. We demonstrated the feasibility of 
this fitting ansatz in CP preserving 2 and 3 generation models using a perturbation 
around the dominant 2-3 sector of the fermion mass matrices. Furthermore, using a 
"downhill simplex " non linear fitting procedure the authors of [18] were able to show 
that even with CP preserving (i.e real) Yukawa couplings the CKM CP violation as 
well as all the other known fermion mass and mixing data could be accurately fit with 
Yukawa couplings which followed the '126 suppressed' dominance pattern. Although 
the Yukawa coupling of the 10 plet Higgs (weakly) dominates that of the 126 in the 
solution of [18], the same is not true of the 120. This is rather different from the 
situation in the perturbative fit obtained by us [IT] where the couplings of the 126 
are very small. Since these two solutions are only early attempts to view the vast 
parameter space of the fit we shall not belabour the difference. 

The authors of [18] fixed 3 of the Yukawa parameters using an ad hoc Z2 symmetry. 
Thus these 3 parameters remain "in reserve" and can be called to fit the values of the 
as yet unknown Leptonic CP phases and Neutrino mass scale besides the actual value 
of 6± 3 (as opposed to a limit). The indications are thus that even a NMSGUT with 
only real Yukawas i.e with a total of 24(= 12 fermion Yukawas + 12 AM Yukawas) 
real parameters (further reduced to 23 by the fine tuning condition) will still be able 
to fit all the fermion data ( 12 masses + CKM phase + 3 CKM angles + 3 PMNS 
angles + 3 PMNS phases = 22 parameters) ! Such a theory will have less parameters 
than even the (unsuccessful) old MSGUT ! Moreover this fitting problem which will 
have 22 data values fit by at most 23 parameters (once the missing neutrino data is 
available and counting phases as "full" parameters !) holds out the enticing prospect 
of predictions regarding 'smoking gun ' baryon violation signals of grand unification. 
Even before such data is available however, it would clearly be preferable to analyze 
the simpler case of spontaneous CP violation first. This shall be a guiding theme in 
our discussions of the NMSGUT parameter space. 

It is clear that that the 120 does not contain any SM singlets and hence the 
analysis of the GUT scale symmetry breaking in MSGUT carries over unchanged 
to the NMSGUT. In particular there is still only one complex parameter(a;) whose 
variation directly affects the vevs and thus the masses in the theory. The additional 
kinetic terms are given by covariantizing in the standard way the global SO (10) 
invariant D-terms [-mO* jk O ijk ] D 
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3 AM Chiral masses via PS 

As in the case of the MSGUT j3l |5] we open up the maze of NMSGUT interactions 
by decomposing 50(10) invariants in the superpotential first into Pati-Salam invari- 
ants and then, after substituting the GUT scale vevs in PS notation we obtain the 
superpotential in the MSSM vacuum in terms of MSSM invariants. The results for 
the old MSGUT case were already given in j3j [5] thus we list only the effect of the 
additional terms in the superpotential. The PS form of Wnmsg is (we have inserted 
line numbers for easy reference) : 



^HiOjki^ijki = k[±i(H^O^ + H^O^ x n (18) 

+ ^{H^o^f Xa& -H^or a ^ {s) ) ( 19 ) 

+H$(d$ m • + Og (z) • (20) 

-\K)° a& K a L (21) 

+^(OW<g+OJ5$JLi) (22) 

-^=H aa (Of a ^; ap + OjfQZp) (23) 

2^2 ^ va u VW s i»' & ' ^ > 

+H a& O a& $] (25) 

^O ijk O ijk = ^[60^0^ + 60^0/^ (26) 

+3(0 (a)(fl) • 0^ {R) + (a)(i) • 0« L) ) (27) 

-60 a "O aA ] (28) 

%O ijm O klm <S> m = |[8z0g0g$/ + 8z0r d 0^ d ^ (29) 

-4v^(0g0 (a)(JI) ■ ^ + 0g0 (8)(L) ■ <&„J L) (30) 

-0^6%, ■ *fa - o$p$ w ■ ( 31 ) 

+4^(0^0^^ - 0°*Ojf $ ; Q/3 ) (32) 

+8(0/^02$^ - 0/^Og$g ad ) (33) 

- 2 (^°- *m*M + Of^Oj^J (34) 
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-2(0 (a)(L) • Og £) - {a)(R) • OW (JO )$ (36) 

-^(O^Og^Wj - Or*0&**<$) (37) 

-2^(0^0^^$^ + O^^f ) (38) 

+4^(0,-0^$^ + Or^Otf*^) (39) 

-2(0//0^ ) $^ ) + OjfdrV*^ (40) 

~^Ad U ^, a p %S {s) ~ U Xa U lxvdt ^(s) ) l 4i J 

-4^(0^^ • (d£ )R X $ A fo) (42) 

+<%) • (<%l x (43) 

JLj[-6v^<(02e{5^ - Og}sW$ A ") (44) 

-12i(OgEr^^2 + OgE^c&j&J (45) 

-12v^(0Wfg w • S„fo - 0gE2 (L) • $^ L) ) (46) 

+6v / 2.(OgE/ aa $^ ) - 0^d<H (47) 

-Bv^C^Eft $g Qd + O a ^EW$£ w ) (48) 

+eV2(o^l P ^ + o^^p^) (so) 

+mor^; aa ^ (si) 

+12i(0 (a)(L) • Eg (£) $^ + W W . Eg w $ A ") (52) 

-^(^^^W - O^Kf^xt) (53) 

-12(OSl)^L) - ^U^W^ (55) 

+6(0f E< V£7 a - Of 2^$^) (56) 

-6v^(0f Er**^ + Ofcr**^) (57) 

+12(0/^E/ Qd ^ - OrE^f/) (58) 

+12(0/-Eg Q/3 $^ - Or^^H) (59) 

-12(Or d ^ )$ (a)cd) (60) 

+«(^ E ^*2f + d^^f x i (6i) 
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+12v / 2[0 (a)(i) • (<&„fo x E« (£) ) (63) 

• x (64) 

+6^(0^2//$^ + O^E.f (65) 

^^(Ogl^ - OgE^/) (66) 

-12^(0«Sr Q $^ + OgE^SL) (67) 

-12v^(OgE (L) • ^l) " Ow E maW • W ( 68 ) 

-6v / 2^(OgE A ^$W^ - Ogs^l-^) (69) 

+6v / 2(O^E( a >(l a + O^SW^W } (7Q) 

-i2v^(or d £, A >;^ + ors,^^) (7i) 

-6V2(O^S^|.^ + A f ^« *#«*) (72) 

-12*0/ a %L$ (73) 

+m(5^ a) ■ + d w L) . ggW^ (74) 

+6i(0^ ) E^$^ ) - CT^t.^^ (75) 

-C Vd/3^Aa («) jxvo.fi \o. (a) J U D J 

-12(0W (L) • $ A ? L) - • * A fo)I$ (77) 

- 6 (°^ S il^( s r - O/^S^w) (78) 

-e^o/sr^ + oftr**^) (79) 

+12(0/"%^ - 0/"%^$/) (80) 

+12(0/Q d S (^ ^A(.)/i _ O^E^W^J (81) 

-I2(0r a ^^ )aa ) (82) 

+ 6 ^(Of; )a/3 S^$gf + Of^E^gf (83) 

+ O^Z?** " K f) (84) 

+12v / 2[0 (a)(R) • ($^ R) x E MW ) (85) 

-^Sl) ' (*A(L) >< ( 86 ) 

-e^gT'S^oWj + &$*tg O^)] (87) 
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The purely chiral superheavy supermultiplet masses can be determined from these 
expressions simply by substituting in the AM Higgs vevs and breaking up the contri- 
butions according to MSSM labels. 

It is again easiest to keep track of Chiral fermion masses since all others follow 
using supersymmtery and the organization provided by the gauge super Higgs effect. 

There are three types of mass terms involving fermions from chiral supermultiplets 
in such models : 



• Unmixed Chiral 

• Mixed pure chiral 

• Mixed chiral-gauge. We briefly discuss the notable features of the mass spec- 
trum calculation and give the actual mass formulae in the Appendix. 



3.1 Unmixed Chiral 

A pair of Chiral fermions transforming as SU (3) x SU (2)l x U(1)y conjugates pairs 
up to form a massive Dirac fermion . For example for the properly normalized fields 

AH,l,-4) = ^± A[l,l,4] = 3|) ( 88 ) 

one obtains the mass term 



2(M + 7](p + 3a + 6u))AA = m A AA 

The physical Dirac fermion mass is then |tba| since the phase can be absorbed by a 
field redefinition. By supersymmetry this mass is shared by a pair of complex scalar 
fields with the same quantum numbers. 

In the MSGUT case there are 19 pairs of chiral multiplets which form Dirac 
supermultiplets pairwise and two Majorana singletons, none of which mix with others 
of their ilk. In the NMSGUT 6 of these pairs become mixed with others of the same 
type leaving 11 Dirac supermultiplets and 2 Majorana supermultiplets (S,Q) which 
are unmixed. If the representation is real rather than complex one obtains an extra 
factor of 2 in the masses. The relevant representations, field components and masses 
are given in Table 1. 



3.2 Mixed Pure Chiral 

For such multiplets there is no mixing with the massive coset gauginos but there is 
a mixing among several multiplets with the same SM quantum numbers. There were 
only three such multiplet types in the MSGUT (i.e R[8, 1, 0], h[l, 2, ±1], t[3, 1, ±§]) but 
in the NMSGUT, as already mentioned there are an additional 5 mixed pure chiral 
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types namely the C{8, 2, ±1], D[3, 2, ±|], fsT[3, 1, ±§], L[6, 1, ±§], P[3, 3, ±f ]. As for 
the multiplet types which had mixed pure chiral mass terms in the MSGUT, the type 
R{8, 1, 0] acquires no new partners and has an unchanged mass matrix since the 120 
has no such submultiplets. However the other two mixed pure chiral multiplet types 
of the MSGUT do acquire new contributions : 

• [1, 2, -l)(h u h 2 , h 3 , h 4 , h 5 , he) 0[1, 2, l)(h u h 2 , ha, h 4 , h 5 , K) = (H«, Sf )a , 

v (15)a nQ n (15)a N m /rr ^( 15 ) V^ 15 ) 5^ D n( 15 ^ 

Here one gets additional an 2 rows and 2 columns relative to the MSGUT since 
the 120-plet contains two pairs of doublets with MSSM type Higgs doublet 
quantum numbers so that the mass matrix 7i is 6 x 6. To keep one pair of light 
doublets in the low energy effective theory, it is necessary to fine tune one of the 
parameters of the superpotential (e.g Mh) so that DetTi = 0. By extracting 
the null eigenvectors of TVO, and T-ChO one can compute the composition of 
the light doublet pair in terms of the doublet fields in the full SO(10) GUT, 
and, in particular, we can find the proportions of the doublets coming from the 
10, 126, 120 multiplets which couple to the matter sector (see Section 5 and 
Appendix C). This information is crucial for investigating whether a fit of the 
fermion data accomplished by using the generic form of the SO (10) fermion 
mass formulae is compatible with the dictates of the MSGUT. In fact precisely 
such considerations led [HI EE] to the conclusion that the MSGUT is Type I 
Seesaw dominated yet gives too small neutrino masses. 

• [3? 1) fa, t 4 , ?5, ij) 0[3, 1, — §](£i, £2, £3, £4, £5, ^6, £7) = (^ 4 )^fa)5 

S fa)> S (R0)' $ 4( J R+)> 0^ S \ Offlo)) ®(HpA, Sp4( a ), £ A 4(a) , Ep4(R0) , ®t(R-)> °p4(s), O^rq)) 

With the contribution of the 120-plet one gets two additional rows and columns 
and the dimension of [3, 1, ±|] mass matrix T becomes 7x7. These triplets 
and antitriplets participate in baryon violating process since the exchange of 
(ti, t 2 , £4, t e , tr)(B(ii,t2, te, tf) Higgsinos generates d — 5 operators of type QQQL 
and luud. The strength of the operator is controlled by the inverse of the t — t 
mass matrix T. 

3.3 Mixed Chiral-Gauge 

Finally we come to the mixing matrices for the chiral modes that mix with the gauge 
particles as well as among themselves. There is no direct mixing between MSSM 
fields contained in 120-plet with gauge particles. However mixing is present via other 
MSSM submultiplets present in MSGUT Higgs fields which further mix with gauge 
fields. This occurs for all such multiplet types except G[l, 1, 0] and X[3, 2, ±|] which 
are unchanged, while for E[3, 2, ±|], F[l, 1, ±2], J[3, 1, ±|] mass matrices acquire ad- 
ditional rows and columns. Thus 
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• [3 5 2, — gK-Ei, E 2 , E 3 , E 4 , E 5 , E & ) © [3, 2, £2, E 3 , E4, E 5 , Eq) 

= l^4i ' ^41 > 9» W 2> ^2 ' A 2 ' U 4i J 93 ^pcA^pca* VpiaV $ W' ^ a1 ' 

The 6x6 mass matrix £ has the usual super-Higgs structure : complex conju- 
gates of the 5th row and column (omitting the diagonal entry) furnish left and 
right null eigenvectors of the chiral 5x5 submatrix E obtained by omitting the 
fifth row and column. S has non zero determinant although the determinant of 
E vanishes. 

. [1,1, -2]( F U F 2 ,F 3 ,F A )©[1,1,2]( F U F 2 ,F 3 ,F 4 ) = (X um /V2, X(r-), 
O m /V2 ) ®(Z m /V2, \ {R+) , 4 yV2) 

The 4x4 mass matrix T has the usual super-Higgs structure : complex conju- 
gates of the 3th row and column (omitting the diagonal entry) furnish left and 
right null eigenvectors of the chiral 3x3 submatrix F obtained by omitting the 
third row and column. T has non zero determinant although the determinant 
of F vanishes. 

• 1? — |](^i) J2, J3, J4, J5) © [3, 1, s\(Ji, J21 J3, J4, J5) 

The 5x5 mass matrix J also has the super-Higgs structure : complex conjugates 
of the 4th row and column (omitting the diagonal entry) furnish left and right 
null eigenvectors of the chiral 4x4 submatrix J obtained by omitting the fourth 
row and column. J has non zero determinant although the determinant of J 
vanishes. 

This concludes our description of the superheavy mass spectrum of the NMSGUT, 
explicit details are given in the Appendix. 

4 RG Analysis 

In [HIE], for the case of the MSGUT, we exhibited plots of the threshold corrections 
(A<5, A\y, A x ) to a G (M x )^ 1 , Sin 2 0w and log 10 M x versus £. Note that we follow the 
approach of [32] in which M x is taken to be the mass of the lightest gauge multiplet 
which mediates proton decay and not the point where the 3 MSSM gauge couplings 
cross. In [15] we synthesized the three batches of information corresponding to the 
three roots of the cubic equation® by exhibiting contour plots of A G , Aw, A x on 
the x— plane at representative values of the other ('slow') parameters A, 77,7, 7. In 
this section we summarize information on the x values allowed by imposing plausible 
'realistic' constraints on the magnitudes of the threshold corrections to the gauge 
couplings. We pay particular attention to the scenario [23] where M x and with it all 
dangerous triplet masses are pushed above 10 l6 GeV. To implement the consistency 
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requirements that the SO(10) theory remain perturbative after threshold and two 
loop corrections and, conversely, that a G not decrease so much as to to invalidate the 
neglect of one-loop effects in the chiral couplings we impose an upper limit of 10 on the 
change in the magnitude of o>q 1 . Similarly we estimate that uncertainty concerning 
the value of the supersymmetry breaking masses translates to an uncertainty of up to 
10% in sin 2 9w(M s ). Finally we expect that the mass of the lightest gauge mediating 
baryon decay bosons should not be lowered by more than one order of magnitude in 
order to respect the current bounds on d = 6 mediated nucleon decay and not be 
raised by more than 2 orders of magnitude so that it remains well below the Planck 
scale. Thus we demand : 

|A G | = (A^Mx))! < 10 
2 > A x = A(Log 10 M x ) > -1 

\A W \ ee \A(sin 2 6 w (M s ))\ < .02 (89) 

We find threshold corrections |5j [H] 

M' 



M x 



A {th) {Log w M x ) = .0217 + .0167^(56; + 3b' 2 - 8b' 3 )Log w 

M' 

M' 

A^ th \sin 2 e w {M s )) = .00004 - .00024^(46; - 9.6b' 2 + 5.6b' 3 )Log w 

M' 



M x 



M' 

A^ th \a~ G }{M x )) = .1565 + .01832 £(56; + 3b 2 + 12%) Log 10 — (90) 

M' Mx 

Where b\ = 167r 2 ^ are 1-loop j3 function coefficients (/% = hgf) for multiplets with 
mass M' . These corrections, together with the two loop gauge corrections, modify 
the one loop values corresponding to the successful gauge unification of the MSSM, 
see [SJ El [T5] for details. 

The parameter £ = XM/rjm is the only numerical parameter that enters into the 
cubic eqn.® that determines the parameter x in terms of which all the superheavy 
vevs are given. It is thus the most crucial determinant of the mass spectrum . 

The rest of the coupling parameters divide into "diagonal" (A, rj, p) and "non- 
diagonal" (7, 7, (, C, k) couplings with the latter exerting a very minor influence on 
the unification parameters. We have therefore fixed the non diagonal parameters at 
representative values ~ 1 throughout. The dependence of the threshold corrections 
on the "diagonal" couplings is also comparatively mild except when coherent e.g when 
many masses are lowered together leading to «g explosion, LogM x collapse or large 
changes in sin 2 8w{Ms). This happens when we lower these couplings too much. For 
lack of space and cogent reason to be exhaustive we give sample plots to summarize 
the most characteristic variation of the threshold effects over the complex x— plane. 
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A crucial point [14] is that the threshold corrections depend only on ratios of 
masses and are independent of the overall scale parameter which we choose to be m. 
Since M x = 10 m - 25+Ax GeV it follows that 



A AfT Mx 



XGeV 

m \ = 10 16 - 25+A * ; |A| GeV (91) 



g^jA\a + w\ 2 + 2\p + uj\ 2 

It is thus clear that this factor will enter every superheavy mass so that they must 
all rise or fall in tandem with Mx he with Ax- 

In Figs. 3,4,5, for the case of the old MSGUT, we give examples of regions 
of the complex a;— plane which are compatible with the unification constraints ( J89l) 
for characteristic representative values of the diagonal couplings. As either of the 
couplings A, rj is lowered below magnitudes of order unity the region of the x-plane 
where the three unification constraints fl89l) are simultaneously satisfied changes in a 
characteristic way. The twin lobes (Fig. 3) in the first and fourth quadrants rotate 
towards the second and third quadrants respectively (Figs 3,4,5). The same behaviour 
is observed in the case of the NMSGUT but with the twist that the dependence on 
the parameter p can induce a considerable broadening of the allowed region when p 
is in a critical region near p c — .15 (Figs 6,7,8) and the other diagonal couplings are 
~ 1. 

An even more remarkable feature of both these sets of plots is that the allowed 
regions of the x-plane are dominantly those where the unification scale is raised above 
10 17 " 25 (jeV(the darkest shaded parts of Figs. 3-8). Thus the NMSGUT seems to point 
unequivocally and generically towards a resolution of the difficulties with d = 5 baryon 
decay and a too low gauge Landau pole by an across the board elevation of GUT scale 
masses. 

As discussed above, due to the successful fit of fermion masses and CP violation 
using only real superpotential couplings, the fitting problem for the real case is even 
more interesting than the general case and certainly more tractable. If the superpo- 
tential parameters are real so is £ and it follows that the 3 solutions X{ of eqn.Q 
form a conjugate pair accompanied by a real solution or else are all independently 
real. We shall show in section 5 that spontaneous CP violation requires that x be 
complex. Due to the presence of a reflection symmetry that interchanges the complex 
conjugate pair of solutions of eqn.flBJ) it is sufficient to study solutions with positive 
real imaginary part only. The complex conjugate pair of solutions exists only for 
£ > —27.916 and are present as a conjugate pair. Then we may ask for what (real) 
values of £ can we obtain complex x values compatible with the constraints of unifica- 
tion. As already seen the answer depends on the values of A, 77, p. In Fig. 9 we show 
a parametric plot (vs £) of the branch £+(£) of the solution of eqn.(jH])with positive 
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Figure 3: MSGUT all couplings ~ 1 : Regions of the x-plane compatible with 
the unification constraints (1891) are shaded. The darkest regions have 2 > Ax > 
1 (corresponding to Mx > !0 17 ' 25 GeV), the next darkest 1 > Ax > the lightest 
shade > Ax > — 1 and the white regions are disallowed. 



imaginary part. The reflection symmetry in the Re[x] axis makes discussion of just 
the positive imaginary part branch sufficient. 

In Figs. 10,11,12 we exhibit 3D plots showing regions of parameter space where 
viable complex x values can be generated from real values of £. This is a first step in 
the goal of determining which regions of the parameter space could actually support 
the successful generic fits of the fermion data discussed above. It is remarkable 
that the combined requirements of spontaneous CP violation and unification leave 
only rather restricted regions of the diagonal parameter space viable. In particular 
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lm W (1=0 12^ ij=0.95 s 7= 07, 7=06) 



3.5 




-3.5 -1.75 1.75 3.5 

Figure 4: MSGUT with A lowered: Regions of the x-plane compatible with the uni- 
fication constraints ( l89|) are shaded. The darkest regions have 2 > Ax > 1 (cor- 
responding to Mx > 10 17 25 GeV), the next darkest 1 > Ax > the lightest shade 
> Ax > — 1 and the white regions are disallowed. 



small values of the diagonal parameters A or r] are not allowed : as may be seen by 
superposing Fig. 9 on Figs. 4, 5 (which are for the MSGUT but are essentially similar 
to the NMSGUT plots for small A, 77) . The reason is that lowering A, 77 results in 
lower thresholds and explosive growth of the gauge coupling. When both A, r\ are 
~ 1, however, one finds that the solution curve shown in Fig. 9 intersects the allowed 
region for £ values from about £ ~ 10 to as low as £ ~ —27.9 (one can approach even 
closer to £ ~ —27.917 below which all solutions are real by tuning p to its critical 
value) . 
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-3.5 -1.75 1.75 3.5 

Figure 5: MSGUT with 77 lowered : Regions of the x-plane compatible with the uni- 
fication constraints ( |89j) are shaded. Further lowering 77 gives a plot almost identical 
to Fig. 4 where the lowering is due to A. 

4.1 Suppression of Baryon decay 

As discussed in the introduction, a very interesting question is whether there are 
viable regions of the parameter space where the theory is still perturbative (i.e 
o>q 1 (Mx) > 10 , the change in sin 2 9w is acceptable and Mx « Mpi anc k) yet 
the masses of the colour triplet Higginos that mediate proton decay are sufficiently 
large as to remove or mitigate the challenge to GUTs posed by the non-detection of 
Baryon decay. Surprisingly-in view of the tight upper bound on the triplet masses 
in the renormalizable SU(5) theory [2T] the answer is generically in the affirmative. 
The rule of thumb of [21] indicates that Mxripiet > 7x 10 16 CreV is sufficient to sup- 
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Figure 6: NMSGUT critical behaviour : p = .01 i.e much smaller than p c = .15, 
all other couplings ~ 1. Regions of the x-plane compatible with the unification 
constraints fl89|) are shaded and are essentially like the MSGUT case. The darkest 
regions have 2 > Ax > 1 (corresponding to Mx > 10 17 ' 25 GeV), the next darkest 
1 > Ax > the lightest shade > Ax > — 1 and the white regions are disallowed. 



press d — 5 Baryon decay operators to acceptable levels. In our case (see Section 
6) there is not one pair but a plethora of triplets -of three different MSSM types 
(t[3, 1, ±|], P[3, 3, ±|], K[3, 1, ±|])- that can mediate Baryon decay. We can use the 
smallest mass of each of these three types as a surrogate for the the SU(5) Mxripiet- 
However in view of the several multiplets involved it is clear that there will be addi- 
tional mixing angles associated with the diagonalization of triplet mass matrices that 
will enter the baryon decay amplitudes. Hence the lower bound on M trip i et is likely 
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Im[x] 
6 



Q_= 1.0, ?j= 95, j>= .7^ 7 =0 6, £= .7 ^ £ = 85, k= .75, m = 8, p=. 15) 




u 1 1 1 ^ FteM 

-6 -3 3 6 

Figure 7: NMSGUT critical behaviour: p = p c = .15. Regions of the x-plane com- 
patible with the unification constraints ( |89j) are shaded , note the broadening of the 
allowed region with merger of the allowed 'lobes'. 



to be somewhat weaker than quoted (say 10 16 GeV for definiteness) . 

In Fig. 13,14,15 we plot the values of Ax, Ac, Ay^ versus £ with other couplings 
chosen at a representative points in the viable ranges indicated by Figs 10-12 and with 
x = (for spontaneous CP violation). There is a sharp peak in m (the overall 

mass scale in all superheavy masses written in terms of the dimensionless vevs) due 
to a peak in Ax at £ = —5. (which on the branch corresponds to x — 1 

On this complex branch this spike is not due to any special gauge symmetry as can 
be checked from the values of the S , O(10)/Gi 23 coset gaugino (i.e E,X,G,F,J type 
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-3 75 -1.8 1.8 3.75 

Figure 8: NMSGUT critical behaviour : p = .51 i.e larger than p c = .15 : Regions 
of the x-plane compatible with the unification constraints ( )89i) are shaded, note the 
reversion to the two lobe pattern of the MSGUT. 



gauginos ) masses which remain distinct and non zero. Rather it is because a certain 
special multiplet, namely the lowest mass eigenstate in the C[8, 2, ±1] sector becomes 
light as one approaches close to to £ = — 5 ( i.e as M x , «g rise). One cannot approach 
the exact value too closely as «g explodes. This is visible in Fig. 16 where we show 
the Mass of this light state versus £. Note however that the 'crevassse' into which 
Mq 171 is to fall is very narrow indeed and that outside the region £ G (—5.5, —4.5) 
both «g and Mq are comfortingly normal. In fact, except exactly at the singularity 
Mc is always above 10 15 GeV A since it too has a factor of m within it. The quantum 
numbers of this multiplet are reminiscent of the intermediate scale colour octet and 
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Figure 9: Parametric plot of x(£) \Im[x(£) > 0, £ G [—27.917,1000). The terminus 
point near (2.3, 0) corresponds to £ — > —27.917 and that near (1, 0) to £ — > oo 



SU{2) triplet multiplets used in raising the unification mass in [23]. Note however 
that we follow the approach the effective field theory based approach of [32J who 
defines the unification scale Mx simply as the mass of the lightest B violating Gauge 
boson. It would be interesting to reanalyze these questions using other approaches 
to the RG flow. The multiplet C[8, 2, 1±] seems to have an analogous effect in this 
approach to that of the triplet and octet thresholds used in [23]. That is to say 
although its mass is still quite large outside the region £ e (—5.5, —4.5) where «g 
and Mx explode, yet it is already lowered sufficiently to raise Mx sufficiently for our 
purposes, as may be seen from Fig. 13. Ax never drops below about +.2, for £ < 
it is always > 1.2 and goes as high as 5.2 when one approaches £ = — 5. Apart from 
a limited special region around £ — 4, A^ is seen to exceed 1 over the entire £ line. 
One unsatisfactory aspect of this approach is that although consistent with MSSM 
unification at 10 16 - 25 GeV it fails to justify it as a physical gauge threshold since even 
the lightest gauge field is heavier than the point where the couplings in the effective 
theory cross in the one loop theory. 

In Table I we give a list of masses and parameter values for which Mx is raised to 
10 18 GeV, while the gauge coupling of the spontaneously broken SO(10) theory is still 
quite small and we have used a value of £ = —8.5 which is far from the singularity. 
In fact approaching closer pushes Mx above 10 18 Ge^. 
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[y=.7,f=.6, £=.7, £=.85, k=.35, m^.8, p=.95} 



£/50 




Figure 10: Values of A, 77, ^' = .02£ compatible with unification stability constraints 
at p = 0.95. 

The composition of the special C[8, 2, ±1] multiplet is worth mentioning. At the 
singular point it is a mixture of only the modes Ci^, Ci,2 which have their origin in 
the (15,2,2)p5 submultiplets of the 126, 126 multiplets. Thus they couple[5] to the 
matter fields as 

W = %y/2f AB {C 1 d A QB - C 2 u A Q B ) + ... (92) 

It is not clear to us if this rather banal coupling is of any phenomenological 
importance. 

Even if one excludes the region close to £ = —5 in order to maintain perturbativity 
and to keep Mx < 10 18 GeV (indicated by the horizontal line at Ax = 1-75 in Fig. 13) 
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z.oy, u.oy 


5[1,3,0] 


2.64 


t[3, 1,2/3] 


3.47, 2.26, 1.44, 1.03, 0.73, 0.18, 0.12 


C/[3,3,4/3] 


2.18 


\/[l,2,3] 


1.69 


W[6,3,2/3] 


2.05 


X[3,2,5/3] 


0.97,0.97,0.6 


F[6,2,l/3] 


0.7 


Zf8 1 21 


2.37 


Parameterl^a/wes 

Mm[Mass](C/mis o/l0 18 Get;) 
Max[Mass](C/mts ofl0 18 Gev) 
Threshold Corrections 

GaugeMassesiJJnits ofl0 18 Gev) 


A = 1.0, r, = .95, 7 = .7, 7 = -6, C = -7 

( = .85, k = .35, m = .8, p = .95, m/X = 5.28 x 10 m Gev 

0.06(G[1, 1, 0]), 0.12( L[6, 1, 2/3] , t[3, 1, 2/3] ) 

3.86(0[1,3,2]) 

A x = 1.739, A w = -0.149, A G = 0.4994 
a G = 0.038 

m XE = 1.01, rfi\ F = 0.41, m Ac = 0.67, 
m\j = 0.61, m Ax = 0.97 



Table 1: Example of Masses and couplings favouring high Mx, £ = —8.5 on the 
x = branch for real couplings 
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{y= 0.7, ?=0.6, £=.1, £=.85, k=.35, m t) =.8, ri=.95} 



£/50 




Figure 11: Values of A, p, £' = .02£ compatible with unification stability constraints 
at 7] = 0.95. 

it is clear that negative £ alone raises the Unification scale by more than an order 
of magnitude. This implies that in a region near £ = — 5 the complex solutions 
are compatible with perturbative unification yet raise the Unification scale to values 
as close as desired to Mpi anc k. To maintain perturbativity, and the Weinberg angle 
prediction and allow neglect of gravitational radiative corrections it is sufficient to 
demand A x < 2 so that M x < 2 x 10 18 GeU ~ 10 _1 M P . From Fig 14 we see that 
the value of aa rises as one raises Mx but this is true only in a very narrow region 
around £ = 5 which is in any case excluded by too large a value of M X - Otherwise 
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{A=1.0, y=.7,T=.6, £=.1, £=.85, k=.35, ^=.8} 



£75 




Figure 12: Values of p,r),£' = .02£ compatible with unification stability constraints 
at A = 1.00. 

the change in «g from the one loop value is quite minor so that perturbativity in the 
gauge coupling is not an issue until one goes very close to the Planck scale. Obviously 
the d — 6 i.e gauge mediated baryon decay operators operators will be completely 
suppressed when Mx is raised towards the Planck scale (say 10 18 GeV ) which will 
raise the d — 6 operator mediated lifetime by 8 orders of magnitude above the 10 36 
years usually quoted for d = 6 processes in supersymmetric GUTs.In Fig. 15 we 
observe that the change in sin 2 6y/ only restricts one from going too far from the 
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{ A=1.0, Y}= 0.95, y= 0.7, 7=0.6, g=.T, g=.85, k=.35, m<>=.8, f>=.95) 

3 r 
2.5 - 

2 - 



Jf 1.5 - 
1 - 
0.5 - 

- 

-20 -10 10 20 
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Figure 13: Plot of Ax against £ on the CP violating solution branch x+(£) at repre- 
sentative allowed values of the diagonal parameters. 

singularity! 

To see what happens with d = 5 baryon decay operators we show in Figs 17,18,19 
the masses of the lightest Baryon decay triplets of each of the three independent types 
t, K, P (see Section 6). 

The mass of the common B-decay triplet £[3, 1,±|] is seen to always lie above 
10 w GeV and is even higher for most of the £ line. We emphasize that the effective 
mass is greater since there will be a rotation matrix element between the lightest 
eigenstate and the modes that couple to the baryon violating d — 5 operators(see 
Section 6). The situation with the P[3,3,±|] and i^[3,l,±|] type triplets is even 
better since they mostly lie above 10 17 GeV !! 

There is however a price to pay for all this mass scale raising and it is enforced 
by the other arm of the Baryon violation- Lepton violation seesaw or balance that 
operates in SO(10) theories. Namely (Fig. 19) the 126 vev also rises above 6 x 
10 16 GeV for most of the £ line. If one wishes to have right handed neutrinos much 
lighter than the GUT scale this would indicate that solutions of the fermion mass 
fitting problem with very small values of the 126 Yukawa couplings are preferable to 
those that rely on inordinately skewed values of the Higgs fractions. Interestingly the 
semi-realistic 3 generation fits we obtained by perturbing around the 2-3 generation 
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( 1=1.0, r)= 0.95, y= 0.1, 7=0.6, g=.1, £=.85, k=„35, ra =.S, p=.95) 
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Figure 14: Plot of A G against £ on the CP violating solution branch at 
representative allowed values of the diagonal parameters. 

10 — 120 indeed have very small Jab [Tfil [T^ (see Section 5). On the other hand 
the values of / are somewhat larger (~ .1 )in the case of the sample realistic fit 
obtained using a numerical fitting procedure [18]. This may be an argument in favour 
of realistic fits close to our semi-realistic ones (which can be searched for [35] by 
starting the downhill simplex procedure near the solutions we found [TT]). If indeed 
the viable solution finally has high 126 vev accompanied by the very small 126 
Yukawa coupling required by the Type I Seesaw and the 10, 120 dominant charged 
fermion fit, it will be yet another instance of the intricate balance of the Fermion 
mass hierarchy and its ouroborotic link to the structure of the apparently completely 
remote GUT scale symmetry breaking and threshold structure. 

Thus we see that the NMSGUT deals with the difficulty regarding d = 5 medi- 
ated baryon decay in a very natural and unforced way. No fine tuning of Yukawas 
or artificially arranged cancellation or any introduction of a plethora of uncontrol- 
lable non-renormalizable terms is required at all. The suppression is generic in the 
viable regions of the parameter space and is practically inescapable on the complex 
branches of the solution of eqn.fjB]) which are structurally imposed requirements when 
CP violation is spontaneous. 
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( A=1.0, rf= 0.95, j= 0.7, f=0.6, g=.7, £=.85, k=.35, ra =.8, p=.95J 




4.2 Capping the Landau Pole 

Furthermore, raising Mx to values near the Planck scale alleviates an even grosser 
difficulty of the MSGUT one loop unification scenario with 210 © 126 © 126 Higgs 
system without threshold effects [22] : the presence of a Landau pole in the gauge 
coupling evolution at Ax ~ 5Mx i.e lying just above the perturbative unification 
scale but well below the Planck scale. The problem is only worsened by the intro- 
duction of the 120. If the Unification scale is raised close to the Planck scale by the 
threshold effects that we have calculated then the strongly coupled dynamics at Ax 
occurs at or close to the Planck scale itself. This strengthens the heuristic arguments 
[25] that envision a (calculable) UV condensation of coset gauginos in the supersym- 
metric GUT (along with chiral vevs related by the Konishi anomaly and therefore 
an associated induced gravity a la AdS-CFT [28J ) which drives the breaking of the 
GUT symmetry |25j. It is worth remarking that such a scenario naturally overcomes 
all the fundamental objections[27] that led to an abandonment of the induced gravity 
scenarios of the 1980s. Moreover since the cutoff of the perturbative theory is about 
an order of magnitude less than the scale of UV condensation and the Planck scale 
- which coincide, it is natural to surmise that the gauge strong coupling dynamics 
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( 1=1. B, Y]= 0.95, y= 0.1, 7=0.6, £=.1, J=.85, k=.35, ma=.8, p=.95) 
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Figure 16: Plot of M™ n against £ on the CP violating solution branch at 
representative allowed values of the diagonal parameters. 

induces gravity characterized by Mp ~ Ax- In this picture the MSSM Grand Desert 
evolution finds 5*0(10) completion when it crosses the superheavy mass thresholds 
and then 5*0(10) quickly defines its own physical UV cutoff : Ax ~ 5M X - A su- 
per symmetric theory with a physical cutoff escapes the objections (raised on grounds 
of ambiguity of cutoff dependent contributions [2 7J) against gravity induced by gauge 
theory dynamics. The coincidence of the scale of condensation and the Planck scale 
is of course the nub of the matter. Previous attempts to construct induced gravity 
from asymptotically free theories had no plausible reason at all why a large Planck 
scale should be induced by gauge theory(e.g QCD or some variety of Technicolour) 
without any intrinsic truly large scale in it at all. Here however we are 'gifted' with 
coincident Planck and strong supersymmetric (therefore holomorphically controlled 
and calculable [25] ) condensation scales with no extra assumptions at all. Combined 
with the enlistment |23j of the MSSM singlets internal to the GUT in the cause of 
specifying the softly broken supersymmetry low energy theory, the whole scenario 
sounds almost too much like a free lunch to be true ! 
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( 1=1. U, rf= 0.95, y= 0.1, 7=0.6, g=.7, £=.85, k=.35, mo=.S, p=.95) 
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Figure 17: Plot of M™ n , the lightest triplet of MSSM type [3, 1, ±§], against f on the 
CP violating solution branch £+(£) at representative allowed values of the diagonal 
parameters. Units of 10 16 GeV. 



5 Fermion Mass Formulae 

As in the case of the MSGUT one imposes the fine tuning condition DetTi = to 
keep a pair of Higgs doublets H^^Hm (left and right null eigenstates of the mass 
matrix 7i) light. The composition of these null eigenstates in terms of the GUT 
scale doublets then specifies how much the different doublets contribute to the low 
energy EW scale symmetry breaking. In the Dirac mass matrices we can replace 
< hi >^ OtiV u , < hi >— > (XiVd- The fermion Dirac masses may be read off the 
decomposition of 16 • 16 • (10 © 120 © 126) given in [31 [5] and this yields [15] 



m u = v(h + f + g) ; r 1 = — cot/3 ; r 2 = — cot/9 

ax «2 

m v = v(h - 3/ + (r 5 - 3)g) ; r 5 i/X 



«6 + iy/3a^, 

m = v(rxh + r 2 f + r 6 g); r 6 = —= — cot (3 

a 6 + tVoa 5 
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( A=1.D, 1J= 0.95, y= 0.7, 7=0.6, £=.7, £=.85, k=.35, m<j=.8, p=.95) 
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Figure 18: Plot of M^ in , the lightest triplet of MSSM type [3, 3, ±§], against f on the 
CP violating solution branch £+(£) at representative allowed values of the diagonal 
parameters. Units of 10 16 GeV. 



m 



9 



v(nh - 3r 2 f + (r 5 - 3r 6 )g); r 5 



4iv3o!5 



■ cot ft 



(93) 



2ig]J^(a 6 + i\^3a 5 ) sin ft ; h = 2\[2ha\ sin ft ; / = — 4^j^ifa 2 sin ft 



The Yukawa couplings of matter fields with 120 Higgs field give no contribution to the 
Majorana mass matrix of the superheavy neutrinos v A so it remains M AB = 8\/2/ab0". 
Thus the Type I contribution is obtained by eliminating v A 



M 



W 

u(I) 
AB 



-M v AB v A v B + VAm V AB v B + 



-M AB 'v A v B + 



-{{m v ) T {M D Y 1 m v ) 



AB 



(94) 



As shown in[T4"l [T5] it is likely that the Type II seesaw contribution is subdominant 
to the Type I seesaw. However the consistency of the assumption that it is negligible 
must be checked and quantified so we also evaluate the tadpole that gives rise to the 
Type II seesaw since the 120 plet does contribute new terms. 
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( A=1.0, rf= 0.95, y= 0.1, 7=0.6, g=.7, £=.85, k=.35, m<j=.8, p=.95) 
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Figure 19: Plot of M% in , the lightest triplet of MSSM type [3, 1, ±|], against f on the 
CP violating solution branch £+(£) at representative allowed values of the diagonal 
parameters. Units of 10 16 GeV. 



For computing O(10, 3, l)i2g vev > inspection of the mass spectrum (Appendix A) 
yields the relevant terms in the superpotential as 



W$ M = M 0-0-^H a ^^O al3 -^=H a ^O al3 



So 



V2 

— + (— ^=ai + 2iv / 377«2 — Vs(a e + i(a 5 )aAvl (95) 
a/2 
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{ A=1.0, rt= 0.95, y= 0.7, 7=0.6, £=.7, £=.85, k=.35, . 8, p=.95) 




Figure 20: Plot of the vev <r of the 126 field responsible for the righthanded neutrino 
mass, against £ on the CP violating solution branch at representative allowed 

values of the diagonal parameters. Units of 10 16 GeV. 



< CL >= 



—7=ai + 2iV3r]a 2 - V3(a 6 + i(a 5 )a 4 {-^j- 
v 2 -/Wo 

and Mo can be read off from Table I to be Mo = 2(M + r](3a - 
neutrino mass is then simply M V AB = 16ifAB < O- >. 



(96) 

p)). The Type II 



5.1 Spontaneous CP violation and the Higgs Fractions a>i,a>i 

As is apparent from the fermion mass formulae the coefficients oti, on are critical for the 
phenomenology of these models. They are calculated by determining the null left and 
right eigenvectors of 7i (the [1, 2, ±1] mass matrix) and can be used to to verify the 
compatibility of the NMSGUT with the realistic generic fits now achievable [TBI \T7 \ [18] . 
To this end we list explicit expressions for the c^, oti in Appendix C. An immediate 
application is to check the conditions under which fits of the fermion data like the 
spontaneous CP violating generic fit of [18] can be realized in the NMSGUT. Intu- 
itively it is evident that complex GUT scale vevs should require complex values of x. 
One finds that the six independent phases [18] that appear in the generic spontaneous 
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CP violation case in 10 — 120 — 126 are given in terms of our quantities by (in the 
convention where a,\ = a± are both real so that the contributions of the lOplet to 
all Dirac masses are real and the phase of a in the Type I seesaw formula has been 
absorbed by redefining the neutrino fields(as also in [18]) ): 

Cu = Arg[a 2 }-^ ; £ u = Arg[ia 6 - V3a 5 ] 

Q-(u = Arg[—\ ; & - = Arg[ — =— ] 

£l - £ u = Arg[r 7 ] = Arg[ 3a 6 + n/3«5 j . ^ _ ^ 

From Appendix C one can check that for real superpotential parameters and real 
values of x the a^ot^ i = 1...5 are real while ae,aQ are pure imaginary. Then it 
immediately follows that except the trivial 126 phase convention dependent phase 
values ( u = Q = — f (essentially from the factor of i that accompanies Jab within 
the parameter Jab) & H other phases are zero and there is no CKM CP violation. 
Clearly the only way to get non trivial phases in the CKM matrix while retaining 
real superpotential parameters is for x to be complex. 

The generic fits 10 — 120 — 126 fits to the fermion data are to be carried out using 
eqns. fl93l) and the Type I seesaw masses [TH [TTJ, [18] . They thus yield information only 
on the ratios 
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«6 
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«5 
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= 7*1 tan (3 ; 


OLx 


a 5 


ir 5 


ct 6 


VS(r 5 - 4) 


a 5 


if 5 


a 6 


VS(r 5 - 4) 


a 6 


(f 5 - Ar 6 ) 


Q!6 


(r-5-4) 



a 2 n 

— = r 2 tan p 
a 2 



tan (3 



tan/3 (98) 

The fitting procedure reduces our ignorance concerning the embedding of the 
MSSM in the MSGUT considerably and thus the stage is set by the formulae and 
analysis given here and in[T5l [TBI [T7l [T8~] for better informed estimates of the exotic 
processes such as Baryon violation and Lepton flavour violation. The latter may 
be crucial for falsifiability tests if the high-Mx route to suppress B-decay(and solve 
several other core problems !) is adopted. The critical unknowns are now the Leptonic 
CP violating phases but the theory is already sufficiently well specified as to permit 
a survey of the possibilities [35] . 
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5.2 Sample fits 



To illustrate the kind of realistic fit obtainable using the 10 — 120 — 126 Higgs system 
we translate the tan {3 = 10 fit obtained in [18] into our notation and also give our 
own semi-realistic fits [T7J for comparison. These expressions also allow us to easily 
see the special features required to ensure the consistency of the fit in terms of their 
implications for the a,, a^. We emphasize however that these fits are significant purely 
as proof of principle. It may also be that the fits with much smaller values of 126 
Yukawas as found by us [17] may be preferred phenomenologically since they allow 
lowering of the righthanded neutrino mass even when the 126 vev has been raised 
along with Mx- Relatively low right handed neutrino masses are often required 
required by Leptogenesis scenarios. This possible dovetailing of the requirements of 
baryon decay(suppression) and (left and right) neutrino phenomenology is perhaps 
another tantalizing glimpse of the ouroborotic balance [16] visible in the construction 
of MSGUTs. 



5.3 Grimus and Kiihbeck realistic fit 



f = e 



-0.0517654i 



,3.95591i 



Dia^.000375, -.021266, .58331} 

/ -0.00016545 0.00727765 

0.0210717 

^ 0.00727765 -0.105497 

/ .00031077 

-.00031077 0.0015122 

V -0.0015122 

r x = .0109799 ; r 2 = 0.00309973 + 0.00129256i 
r 5 = 291.223 - 306.971? ; f 5 = 0.19026 - 0.134289i 
r 6 = 0.110976 - 0.0852565i (99) 

From these one easily deduces 



ai - n tan/3 = 0.1098 ; — = r 2 tan/? = 0.03358e- 395i 



«1 «2 

e l-578i 5gll2 ( 10Q ) 



a 5 Zr 5 1.578i 



«6 VS(r 5 - 4) 

= e ir ' 48l .003198 



nr. zr 5 tan/3 , - 7JS . 



"6 y/3(r 5 - 4) 
«. = (f 5 -4r 6 )tan/3 =e _ 1 ,3 6 ^ 00778 
a 6 (r 5 - 4) 
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In [18] the vev (wr ) of the 126 is taken somewhat arbitrarily to coincide with the 
GUT scale 2 x 10 16 GeV. This could lead to inconsistencies since there is absolutely 
no reason why wr can be assumed to take the fixed value 2 x 10 16 GeV '. 

In [18] the increase in the size of neutrino masses is traced to the fact that f^ 1 
is large and that a certain coefficient has a large magnitude : \r^\ = 3008.9 ~ 
|(r 5 — 3)/r 6 |. Since the denominator contains the two small numbers 26., sl we see 
that the largeness of rp arises from the requirement that the fraction of the MSSM 
doublet H contributed by the 120 plet is very small. The same is evidently also true 
for the 126 since — is also very small. It remains to be seen if these requirements 
are all in fact realizable in the NMSGUT or whether a different type of fit such as 
the one with much smaller 126 Yukawas which leads to very different values of the 
ratios of a's is more successful. 

The fit of [18] was accomplished by fixing f 12 = f 2 z = Qiz = by imposing a Z 2 
symmetry and also by assuming the vev of the 126 to be Mx = 2 x 10 16 GeV. Since 
this vev enters only in the overall size of the neutrino masses it is clear that some 
leeway regarding that overall scale is still present if this assumption is not made given 
that they used the remaining parameters to fit the measured mass squared splittings 
of the neutrinos rather than the masses themselves. Moreover in the fit of [18] the 
leptonic CP phases emerge as (forced) outputs and have the values : 



8 P = .21 ; of = -3.1239 ; a 2 = -3.1404 (101) 

Since these phases have not yet been measured, and there no reason to assume 
that the Z 2 symmetry and the assumption regarding the 126 vev have any basis in 
fact, it follows that we should regard these 4 parameters i.e / 12 , fa, gi3, ojMx as 
lying in reserve for the day when the Leptonic CP angles and the neutrino mass scale 
are finally measured. This equality of the number of remaining unknowns and the 
number of remaining free parameters marks a rather poignant conjuncture for the 
NMSGUT leaving it poised on a razor's edge between inconsistency with the data 
and a grand leap towards an actual prediction of the rate of exotic processes such as 
baryon decay based firmly on an extrapolation of measured data and clear theoretical 
assumptions. 



5.4 Semi-Realistic fits with ultra small 126 Yukawas 



r 5 



.016372 ; f 2 = 0.98039 ; r 6 



2.67 



f 5 = 0.079405 



0.0187 



/ -0.0002988 
0.0000394 
\ -0.0222691 



0.0000394 -0.0222691 
0.009291 -0.134823 
-0.134823 0.481333 
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./ 



/ 0.0000132 
0.0000277 
-0.0001 





-0.00164 
0.02478 



0.0000277 -0.0001 
-0.0000146 0.000057 
0.000057 0.000040 

0.001644 -0.02478 
0. -0.119710 
0.119710 



(102) 



(103) 



The eigenvalues of h, f,g are 



h 
f 
9 



.518 ; 
1.33 x 10 
±0.122 ; 



.0277 



1.92 x 10 



-5 



-1 



1.11 x 10 



-4 



0.171 x 10' 







(104) 



Observe that the values in f AB are much smaller than the fit of [18]. It is also 
easy to check that the associated aij, 6ti values are very different. It remains to be 
seen which type -if either- is actually realizable in the fully specified NMSGUT that 
we have described in this paper. 



6 d = 5 Operators for B, L violation 

In [5] we worked out the effective d = 4 superpotential for B + L violating processes 
due to exchange of colour triplet superheavy chiral supermultiplets contained in the 
10, 126 Higgs multiplets. These included a novel channel due to decays mediated by 
exchange of triplets tu\ contained in the 126 Higgs irrep. Evidently the inclusion of 
120plet Higgs will lead to additional channels for baryon violation. These can be eas- 
ily derived using the Pati-Salam decomposition of the 16.16.120 SO(10) invariants[3j 



W° M = 2V2g AB [h 5 (d A Q B + e A L B ) - h 5 (u A Q B + v A L B )\ 

-2V2g AB [V2L 2 ■ (QaQb)' + F A L A L B + y/2t e Q A L B 
+2-y/2L 2 • [u A d B )' + y/2t 6 (u A e B - d A i> B ) + 2F 4 £> A e B ] 
-2V2g AB [2C 3 2 A Q B - 2C 3 {i A Q B + -j=h 6 (d A Q B - 3£ A L B ) (105) 
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--^h(u A Q B - 3u A L B ) + 2D 3 e A Q B - 2E§v A Q B 

+2E 6 d A L B - 2D 3 u A L B ] - 2iV2g AB [eJ 5 d A d B 
+2K 2 d A e B - eK 2 u A u B - 2J 5 u A u B 

-V2ei 7 d A u B - \p2t-j{d A v B - e A u B )} - 2g AB [eP 2 Q A Q B + 2P 2 Q A L B ] 

We have suppressed G 32 i indices and used a sub-multiplet naming convention specified 
in Section 2, the primes in the terms containing colour sextets (L, L) are explained 
in the caption to Table 1. The carets on the fields indicate that they are fields of the 
NMSGUT and are unitarily related to those of the MSSM(M X ). 

In order that the exchange of a Higgsino that couples to matter with a given 
B + L lead to a B + L violating d = 5 operator in the effective theory at sub GUT 
energies it is necessary that it have a nonzero contraction with a conjugate (MSSM) 
representation Higgsino that couples to a matter chiral bilinear with a B + L different 
from the conjugate of the first B + L value. On inspection one finds that not only 
the familiar triplet types [3, 1, ±|] C 120 i.e {t( 6 ), t( 7 )}[3, 1, |] and {£(6), £(7)} but also 
the novel exchange modes from the P[3, 3, ±|] and K [3, 1, ±|] multiplet types can 
contribute to baryon violation. In the case of the 126 the P\,Ki C 126 multiplets 
did couple to the fermions but P\,Ki C 126 did not. The 120 however contains 
both P 2 , P 2 and K 2 , K 2 . Since these mix with Pi, P\ and K\,Ki, a number of fresh 
contributions appear. 

The multiplets P 2 [3, 3, — |], P 2 [3, 3, |], K 2 [3, 1, — |], K 2 [3, 1, |] satisfy the require- 
ment regarding B + L quantum numbers of the fields they couple to. Note in par- 
ticular that these novel exchanges always lead to contributions in which at least one 
and possibly both pairs of final state family indices are antisymmetrized. 

On integrating out the heavy triplet Higgs supermultiplets one obtains the follow- 
ing additional effective d = 4 Superpotential for Baryon Number violating processes 
in the NMSGUT to leading order in mw/Mx- We have taken the opportunity to 
insert a missing overall sign and correct minor sub/super-script typos in [5] : 



W £ff^° = - l abcd(^Q a Q b Q c Ld) - R AB cD^e A u B u c d D ) (106) 
where the coefficients are 



L A bcd — S^hj^shcD + S^h AB fcD + S 2 f AB hc B + S 2 2 f AB fcD 

— Sih AB gcD — S 2 f AB gcD + ^f2^P^ X ) 2 9AciBD 

- {V- v ) 2 2 ~g AC ~g BD (107) 

and 

Rabcd = S^h AB hcD — S 2 h AB j C D — S 2 f AB h C o + S 2 IabIcd 
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cd + iV%S 4 2 JabJcd 

+ $6 9ABh>CD — iS-^gABhcD ~ S§g AB f CD + ^7 2 9AB f CD 

+ iS-JJiABgcD - iS 2 7 f AB g C D + V2S 4 7 f AB g CD 

+ iS 6 7 g AB g C D + S 7 7 g AB g CD - \Pl{Ki~ x )xf ad~9bc 

- {1C- l ) 2 2 ~g A D~gBc (108) 

here S = T _1 and T is the mass matrix for [3, 1, ±2/3]-sector triplets : W = 
PT/tj + while 

h AB = 2V2h AB f ab = iV2f AB 9ab = ^9ab (109) 

The fields in terms of which we actually require the decay operators are those of 
the MSSM which do not carry carets and are unitarily related to the GUT fields[l4J: 

Qa = Q* A ' A Qa' ; L A = C* A t A L A i ; u A = {V u )* A , A u A > 

~d A ={V d y A , A d A , ■ I A = V e \, A e A , (110) 

The Unitary MSSM C MSGUT embdedding matrices Q,L,V u ' d ' e are fixed 
partly by the unitary ambiguity in the GUT fields which allows one to change 
h,f,g — > U T (h, f, g)U at will, and partly by fitting the fermion mass data to the 
corresponding GUT expressions [141 [T5l [T6l ITT] . This fitting procedure determines 
the embedding matrices up to certain phases and furnishes expressions for the GUT 
couplings in terms of the fermion mass data of the MSSM(Mx) and the embedding 
matrices. So far solutions of the highly non-linear fitting problem have been found ei- 
ther by expanding around the dominant 2-3 sector of the fermion hierarchy |15 | [TBI ITT] 
or by making ad-hoc simplifying assumptions concerning the Yukawa textures to re- 
duce the number of unknowns and then using a non-linear fitting method such as the 
"downhill simplex" [18J. The complicating factor is that we still lack some important 
information on the overall mass scale of neutrinos, leptonic mixing (0^ 3 ) and leptonic 
CP phase information 5 P , af, af. Till they are available one can at most survey the 
possibilities[35j 

7 Discussion and Outlook 

In this paper, motivated by successful fits of the fermion data[TTl[T8] which evade the 
difficulties that forced an abandonment [T4"l [T5] of the hope[9] that the 10 , 126 FM 
Higgs system would be sufficient to describe the entire fermion mass spectrum, we 
specified the ingredients of a New Minimal Supersymmteric GUT based on the gauge 
group 50(10) and the 210 © 10 © 120 © 126 © 126 Higgs System. While inherit- 
ing the Higgs system responsible for GUT scale symmetry breaking unchanged from 
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the MSGUTpU El IE], and fitting all the fermion data successfully, the NMSGUT 
is able to describe CKM CP violation [T8] even when it carries the handicap of real 
superpotential couplings. Since this handicap reduces the number of free parame- 
ters drastically it results in a theory with one less free superpotential parameter than 
even the original MSGUT! Moreover several parameters are yet available to accommo- 
date the new information associated with the as yet unknown Leptonic CP violation 
phases, 9^ angle and the over all neutrino mass scale. If this information became 
available it would determine the unknown but crucial embedding matrices describing 
the relation between the effective MSSM in the grand desert and the NMSGUT. 

Using the techniques we developed for the MSGUTjU [5] we computed the super- 
heavy spectrum for the NMSGUT and used it to compute threshold effects in the 
gauge evolution. Focussing on the symmetry breaking branch which allows sponta- 
neous CP violation we found that the Unification scale defined as the mass of the 
Baryon number violating gauge fields was always raised above the one loop values 
on this branch. This increase could take Mx to values as large as 10 18 GeV while 
still remaining in the perturbative domain. Thus gauge mediated Baryon decay can 
easily become unmeasurably small in this theory. Together with Mx all other masses, 
in particular those of the three triplet types that mediate d = 5 baryon decay also 
rise and can be taken(effectively) well above 10 16 GeV. Thus not only d = 6 but also 
d — 5 baryon decay- which has persistently been a hobgoblin in the fairy story of 
Grand Unification- can be banished to the periphery -if not completely outside -of 
the measurable realm. 

Increase in Mx to ~ 10 18 GeV required by casting the GUT field in the role of 
supersymmetry breaker ( sugrone ?) calls for a reanalysis of the spontaneous breaking 
and spectral determination using the full Sugry potential. Since the GUT singlets are 
well poised to take over the role normally played by the Polonyi singlet or the hidden 
sector [23], we expect the minimum of the global Susy theory to morph into a more 
rigorous solution which also incorporates supersymmetry breaking and suppresses 
the cosmological constant. The participation of singlets from the 126 which couples 
directly to the matter multiplets implies that there may be some (small) deviation 
from universality in the soft supersymmetry breaking masses that arise. 

The increase of Mx provides resolution of another nagging and persistent difficulty [22] 
in the MSGUT : the Landau pole in the gauge coupling evolution above Mx- Since 
M x is closer to the Planck scale the presence of the SO(10) Landau pole next to 
the Planck scale only strengthens our speculation that the UV condensation to be 
expected in such a supersymmetric Asymptotically Strong (AS) theory [22], [25] acts as 
both a physical cutoff for the perturbative SO(10) theory and the scale of an induced 
gravity that arises from this theory. Much remains to be done to implement this sce- 
nario but we made a beginning in [25] by demonstrating using Supersymmetric strong 
coupling heuristics [26] that in a toy ASSGUT the condensation actually takes place 
and breaks the (toy) GUT symmetry and that the vevs responsible are calculable. 
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It is encouraging that the development of the theory in regard to apparently unre- 
lated features has naturally brought us to the point where a number of intractable 
fundamental features have suddenly become pliable to a synthetic interpretation. 

We also gave complete formulae for the fermion masses and baryon violating 
effective superpotential in the NMSGUT, including very lengthy analytic expressions 
for the Higgs fractions (ctj, oti) which are determined by the GUT parameters(after a 
fine tuning) and are crucial ingredients of both the masses and the d = 5 B- violation. 
We illustrated some aspects of their use when examining the realizability of a generic 
fit within the NMSGUT with real parameters by taking the proof of principle fit of [18] 
as as an example. However we note again that much remains to be done to properly 
map out the space of fermion fits in which the fit of [H] represents only an isolated 
point of no particular sanctity. Indeed, if the semi-realistic fits we found in [17] using 
a quite different method based on an expansion around the 2-3 generation core of 
the fermion hierarchy extend to the realistic case with CP violation, we expect the 
new fits to be distinct from those found in [18] and have much smaller values of 126 
couplings. This distinction would then have important phenomenological implications 
since it would lead to much lower right handed neutrino masses ; as indeed is required 
in many Leptogenisis scenarios. 

In consequence of the lowered rates for d = 5 processes which long represented 
the main hope for detection of Baryon decay and for falsifying Susy GUT models, 
it may not be possible to detect such decays, at least in the foreseeable future. In 
that case Lepton Flavour violation[34|will remain the only arena where this theory 
has predictions that may be testable in the near term. 

In sum, the NMSGUT having inherited the strengths of its parent is revealing new 
virtues that illuminate what was mysterious in the MSGUT and promise to carry the 
slow caravan of Grand unification not only across the Grand Desert that was its first 
horizon but across threshold jungles beyond that first horizon up into the rarefied 
heights where gauge forces and gravity meld into their primordial pleromal|25j unity. 
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Fie\d[SU(3),SU(2),Y] 



PS Fields 



Mass 



A[l,l,4],A[l,l,-4] 

M[6,1,|],M[(6,1,-|] 

7V[6,1,-|],7V[(6,1,|] 

0[l,3,-2],0[(l,3,+2] 



W[6,3,f],W[(6,3,-§] 

J[3,l,f],/[(3,l,-f] 

5(1,3,0] 

Q[8,3,0] 

C/[3,3,|],C7[3,3,-|] 



V[1,2,-3],V[1,2,3] 
S[6,2,|],S[(6,2,-|] 
F[6,2,-|],F[(6,2,i] 
Z[8,l,2],Z[8,l,-2] 



V 44 ^ 

(-R+) ^44(fl-) 



/v=r (R+) v^'A^ \ 

/yr'(-R-) \ 

) Zj (i?+)JA<i ? 

_ ^.44 

S 44(L) 



S PP(L)' S (L) 
^ i?(R+)> 4>A{R-) 

4 



0/x(L)> 04(L) 



£ 44 



(Appal* 



2(M + r](p + 3a + 6cu)) 
2(M + rj(p - a + 2a;)) 
2(M + r](p-a - 2a;)) 

2(M + r?(3a -p)) 

2(M- 77(0 + 3?)) 
-2(m + A(p + a + 4a>)) 
2(m + A(2a -p)) 
2(m - A(a + p)) 
-2(m- A(p-a)) 

2(m + 3A(a + a>)) 
-2(m + A(a> - a)) 
2{m- A(a + a;)) 
2(m + A(p- a)) 



Table 2: (i) Masses of the unmixed states in terms of the superheavy vevs . The 
SU(2)l contraction order is always F a F a . The absolute value of the expressions in 
the column "Mass" is understood. For sextets of SU(3) the 6 unit norm fields are 
denoted by a prime = S^ p , p, > 9, = y/2 and similarly for 6. 



Appendix A : Tables of masses and mixings 

Here mixing matrix rows are labelled by barred irreps and columns by unbarred, 
(i) The masses of 13 Unmixed cases are given as Table II. 

ii) Mixed states 



a) [8,2,-l]( C U C 2 ,C 3 )0[8,2,1]( C U C 2 ,C 3 ) = E**, Of )© 

(KvKvO A ai )(A=l 8) 

/ 2(-M + rj(a + u)) -i(u-p)(\ 

2(-M + r](a + u)) + p)( 

V i(u-p)C i(uj + p)( -m + fa / 
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b) [p, -IK D U D 2 ,D 3 )©[3,2,f]( D,,D 2 ,D 3 ) = (££, Eg, )© 

I 2{M + r](a + u)) (iuj + ip - 2ia)( \ 

2(M + 7]{a + 3w)) (-3icu - - 2ia)C 

\ — (ia; + — 2ia)( (3iu + ip + 2ia)( m + |(a + 2c<j) J 



c ) [3, 2, — |](£'i, £"2, -£"3, E4, E 5 , Eq) © [3, 2, |](i?i, £2, £3, -E4, E 5 , E 6 ) 



3 J ±j -^ci, -^4, -^o, -^oy w [<->,.", 3 j 
1 /24a 

'41 ) "41 ) ^(s)2' ^2 ' "2 LI 



- Sff , <& ^ Af Q , Off) © (E^E^, tf« t , ^ Qi , A, ai , o£ ) 



-2(M + r;(a - w)) 














(ia; — 


ip + 2ia)C 





-2(Af + 7j(o - 3w)) 




2ir\a 


igV2W* 


(— 3iaj 


+ ip + 2ia)( 





2iy/2rj5 


-2(m + A(a - w)) 


-2\/2\uj 


2g(a* - u>*) 









— 2irfa 


-2v / 2Aaj 


— 2(m — Xlu) 


V2g(u>* -p*) 









-igV2~v* 


2g(a* - u;*) 


gV2{oj*-p*) 










(— + ip — 2ia)£ 


(3iaj — ip — 2ia)£ 









-(m - 


L §a - fpw) 



d) [1,1, -2]( F U F 2 ,F 3 ,F 4 )©[1,1,2]( F ± ,F 2 ,F 3 ,F 4 ) = , A (i? _), ^ 

^ E (flQ) flj(15) x 44 \ 

A («+)> 72) 

/ 2(M + 7](p + 3a)) -2iV3rja -gV2a* -6i(u \ 

2iV3rja 2{m + A(p + 2a)) y/2tegu* v^C^ 

-^v 7 ^* -y/2&igw* 

\ 6«Cu; v^Co" m + ap J 



e) [1, 1, 0](Gi, G 2 , G 3 , G 4 , G 5 , G 6 ) ee (0, 0( 15 ), «/>g, -^i, ^±1, ^W^^ ) 



m 


V2 
—irjcr 





m + 2Aa 
2 v / 2Aa; 



V 







—vrjay 2 




V^Acj 

2>/2Aa; 
m + A(p + 2a) 
—irjy/Sa 





^^7a• 
V2 

—ir]^/3a 


M + 77 (p + 3a - Quj) 



—ir\cs 

~7T 

M + r](p + 3a - Quj) 


2 



f ) [1, 2, -1] (hi, h 2 , h 3 , h 4 , h 5 , h e ) ©[1,2, 1] (h, h 2 , h 3 , h 4 , h 5 , he) 



^2 '73"' ^2 '^2 J^l-^aU^ai' 



y,(15) v (15) n n (15)> 



al 



/ 




~yV3(u) — a) 


-7V / 3(w + a) 


—70" 


kp 


—\fZikij) 


\ 




-7\/3(a; + a) 





-(2M + 4»?(a + w)) 







i(p + 2u))C 






7v / 3(^ — a) 


-(2M + 4j/(a-w)) 





~2»?ctv / 3 


V3Cw 


-i(p - 2uj)( 






— 0-7 


-2»7o-v / 3 





— 2m + — a) 




VSi^a 






pk 








— m 






V 


V3iku 


i(p - 2a;)C 


-i(p + 2uj)( 


-V3iCa 




—mo — 


/ 







V5ga-* 
2 

V5ga* 
2 





48 



The above matrix is to be diagonalized after imposing the fine tuning condition 
DetH = to keep one pair of doublets light. 

g) [3, 1, -|](Ji, J 2 , J 3 , J 4 , J 5 ) © [3, 1, |](Ji, J 2 , J 3 , ^4, J 5 ) 



/ 2(M + r](a + p - 2u)) -2rp 2y/2rfa -ig^2a* 2((a - 2u) \ 

2r)a -2(m + \a) -2y/2\u -2igV2a* -a( 

-2V2rja -2^/2\uj -2{m + A(a + p)) -Aigu* y/2a( 

~ig^2a* 2^2iga* Aigcu* 

V 2C(a-2cu) a( ~V2>C m Q + f(p - 2w) 7 



h) [3, 1,11(^,^2) ©[3,1, -§]( K!,X 2 , ) = (Eg +) ,Og +) )©(E /i4(il _ ) ,0 A4(R _ ) ) 



2(M+_^(a + p + 2a;)) 2((a + 2cu) 
2({a + 2uu) m + f(p + 2c«;) 



i) [6,1,-|](L 1 ,L 2 )©[6,1,|]( L 1; L 2 , ) = 0'^))© 

(^/ZP(s)(RO)> O'fiv(s)) 

( 2(M + r](p-a)) -2i(oj \ 
y 2i(uj to — |a y 



j) [3,3,|](P 1 ,P 2 )©[3,3,-|]( P 1; P 2 , )^(S( L 4 )) Of i 4 ) )©(S /i4(L) ,d, 4(L) ) 

/ 2(M + 77(a-p)) 2oC \ 
V 2a( mo - fp / 



k) [M.OKi^i^^V 



/i(iJO) ^ 

(to — Aa) — v^Acj 



^ — \/2\uj to + A(p — a) 
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1) [3, 1, |](*1> t2,t3,t4,t 5 ,ts, t 7 ) © [3, 1, — §](£l, t2, ^3,^4,^5, ^6, h) = (H^, , 
S (o)' S CRO)' $ 4(-R+)' 0^ s \ OfjiQ)) © (^4, S p4(a ) , S A4 ( a) , S A4 ( R0 ) , $^ (jR _ ) , A4 (s) , O/^itt))) 







7(a +p) 


7(p - a) 


2V2«a;7 


icry 


V2ka 


\/2*fc<^ 


\ 




7(P - a ) 





2M 








v/2aC 


\/2iu>C 




t(p + a ) 


2M 





4v / 2«^'7 


2i?7(T 


-a/2< 


V2iuj( 






-2^/2iu)-y 


— 4-\/2icj?7 





2Af + 2r)p + 2?7a 


-2v^2??o- 


2ia;C 


2C,a 






icT'y 







2v / 2r?cr 


-2m- 2A(a + p- 4w) 


V2io-C 


-V2C,a 






y/2ka 




V2aC 


-2i£u> 




mo + #a 


-fpu 




I 


— \f2ikui 




— \/2iu>C 


2(a 


v/2?C 




mo + §p 


/ 



m) [3,2 J |](X lj X 2j X 3 )e[3,2, 



•3] (^1)^2) -^3) 



/ 2(m + A(a + u;)) 
-2v^Acj 
-2#(a* + 



-2V2Aa; 
2(m + Acj) 



-2g(a* +uj*) 




Appendix B : SU(5) x U(l) Reassembly Crosscheck 



The internal consistency of these spectra and couplings can be verified by consid- 
ering special values of vevs, e.g 

p = a = ±LU (HI) 

where the unbroken symmetry includes SU(5)[4|. Then we find that the MSSM 
labelled mass spectra and couplings given in Appendix A do indeed reassemble into 
SU(5) invariant form. If we insert a = —uj = p in the mass matrices of Appendix 
A we find that, after diagonalizing the mass matrices of the submultiplets that mix, 
the resultant spectra group precisely as indicated by the decompositions below with 
all the subreps of a given SU(5) irrep obtaining the same mass and correct phases to 
permit reassembly. The delicacy of this reassembly is a very non-trivial consistency 
check of our results. 



H = 10 = 5 1 + 5_ 1 

5x = / il (l,2,l) + ti(3, 1,-|) 

5_i = ^(1,2,-1) +ti(3,l,|) 

(112) 
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£ = 126 = l_ 5 (G/ 4 ) + 5_i + 10_ 3 + 15 3 + 45i + 50_i 

5_! = ^3(1,2,-1) +t 3 , 4 (3, 1,|) 

10_ 3 = F 1 (l,l,2) + J 1 (3,l,~) + S 2 (3,2,^) 

T5 3 = 0(1, 3, -2) + £ 1 (^2,-1) + 7V(6,1,1) 

45x = / l3 (i,2,l) + * 3 (3,l,-^) + Pi(3,3,~)+^ 1 (3,l,|) + D 1 (3,2,~) 

+ 1,(6, 1,-^ + ^(8, 2,1) 

50_! = A(l, 1, 4) + i 3i4 (3, 1 J) + L> 2 (3, 2, ? -) + W(6, 3, |) + M(6, 1, - j|) + C 2 (8, 2, -1) 

_ _ (H3) 

£ = 126 = l 5 (G 5 )+5 1 + 103 + 15_3 + 45_i + 50 1 

5i = /i 2 (l,2,l) + t 2 , 4 (3,l,~) 

T0 3 = F ± (l, 1, -2) + Ji(3, l,l) + £ 2 (3,2,-l) 

15-3 = 0(1, 3, 2) + ^(3,2,^ + ^(6, 1,-1) 

45_! = ^ 2 (l,2,-l)+f 2 (3,l,^) + P 1 (3,3,^) + X 1 (3,l,-^) + J D 1 (3,2,l) 

+ ^(6,1,^ + ^(8,2,-1) 

50x = A(l, 1, -4) + i 2>4 (3, 1, ~) + L> 2 (3, 2, -I) + W(6, 3, -^) + M(6, 1, jj) + C 2 (8, 2, 1) 

(114) 

$ = 210 = 1 + 5_ 4 + 5 4 + 10 2 + 10_ 2 + 24 + 40 2 + 40„ 2 + 75 

lo — Gl,2,3 

5_ 4 = fe 4 (l,2,l) + t 5 (3,l,-|) 

5 4 = ^(1,2,-1) +t 5 (3,l,|) 

10 2 = F 2 (l, 1, 2) + J 2)3 (3, 1,-^+^,4(3, 2,1) 

T0_ 2 = F 2 (l, 1, -2) + J 2 , 3 (3, 1,1) +^3,4(3, 2,-1) 

24 = (1, 1, 0)^,2,3 + 5(1, 3, 0) + X 1>2 (3, 2, -^) + X 1>2 (3, 2, jj) + i? 1;2 (8, 1, 0) 

40 2 = V(l, 2, -3) + £ 3 , 4 (3, 2, 1) + J 2 , 3 (3, 1, -1) + C7(3, 3, -1) + Z(8, 1, 2) + Y(6, 2, 1) 
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Due to the 120plet one obtains the additional SU(5) invariant mass terms: 

(mo + pp)5 5 + (m + pp)10 o 10 o + (m - ^p)45 45 

+2kp (5 5 H + 5 5 ff ) - 2v / 3p(C 5 5 s + C5 5 2 ) + 2(Cor5 5* + C^5 5$) 

+6ip(C10 o T0s + CTOolOs) + ^(C^lOoTOci, + Co"T0 o 10$) 

+2p(C45 45 £ + C45 45 E ) (117) 

Where every SU(5) invariant has been normalized so that the individual G123 
sub-rep masses can be read off directly from the coefficient of the invariant for com- 
plex SU(5) representations which pair into Dirac supermultiplets and is 2 times the 
coefficient for the real representations which remain unpaired Majorana/Chiral su- 
permultiplets. 



(115) 



40_ 2 = V(l, 2, 3) + £3,4(3, 2, ~) + J 2 , 3 (3, 1, ^) + 17(3, 3, ^) + Z(8, 1, -2) + Y(6, 2, ~) 

10 - - 10 5 - 5 

75 = (l,l,0)G 1 , 2 ,3 + /(3,l, y ) + /(3,l,- y )+X 1)2 (3,2,--)+X li2 (3,2,-) 

+ 5(6, 2, |) + 5(6, 2, ~) + i2 1>2 (8, 1, 0) + Q(8, 3, 0) 

O = 5i + 5_i + 10_ 3 + 10 3 + 45i + 45_i 

5! = /i 5 ,6(l,2,l)+i 6 , 7 (3,l,--) 

5_i = ^5 >6 (l,2,-l) + t 6 , 7 (3, 1,|) 

10_ 3 = F 4 (l,l,2) + J 5 (3,l,-^) + 5 6 (3,2,i) 

T0 3 = F 4 (l, 1,-2) + J 5 (3, 1,^+^(3, 2, -1) 

45x = ^ 6 (l,2,l)+t 6>7 (3,l,-^) + P 2 (3,3,-^)+X 2 (3,l,^)+ J D 3 (3,2,-^) 

+ L 2 (6,l,-^) + C 3 (8,2,l) 

45_! = ^ 6 (l,2,-l) + f 6J (3,lJ) + P 2 (3,3,^) + X 2 (3,l,-^) + J D 3 (3,2,^) 

+ L 2 (6,lJ) + C 3 (8,2,-l) 



(116) 
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Appendix C : Doublet fraction Coefficients a i: ai 

In this appendix we give the explicit expressions for the coefficients ctj, ojj obtained 
by first imposing the condition DetTi = and then solving the equations to determine 
the normalized left and right eigenvectors of 7i. 



m 9 X 
mo = cr 



m 



N = 



N = 



(1 -3x)x(l + x 2 ) 
\\ (l-2x + x 2 ) 

1 



\/l«il 2 ^ 


- «2 2 "f 


- d 3 2 + d 4 2 4 


- |« 5 | 2 4 


- "6 2 






1 







Y^laip + \a 2 \ 2 + \a 3 \ 2 + \a 4 \ 2 + \a 5 \ 2 + |«6| 



A = N{a 1 , a 2 , a 3 , a 4 , a 5 , a 6 } = {a u a 2 , a 3 , «4, oc 5 , cte] 
A = N{ai, a 2 , a 3 , a 4 , a 5 , a 6 } = {&i, «2, «3, «4, «5, 

= 8^! = (ri 2 r] 2 XP + ( 2 ( 2 XP 1 +ri ((r ] XP 2 + ((r ] XpP 3 + m ri 2 XpP 4 + r] 2 Xp 2 P 5 ) 

&2 = (m 2 7 r ? AQo + ^oC7CAgi + 7^oC 2 AQ 2 + A;CC 2 AQ3 + C 2 7C 2 Ag4 + 7CC 3 AQ5 

+km (r]XQ 6 + ('y(XpQ 7 + 'y( 2 XpQs + rn 'yr]XpQ 9 + k(i]XpQi + l r] X p 2 Q X1 ) 

a 2 = {m o 2 -ir]XQ + m o C,-iC^Qi + 1 m ( 2 XQ 2 + k ( ( 2 XQ 3 + ( 2 ( 2 XQ 4 + ^ ( ( 3 XQ 5 

+k m ( rj X Q 6 + C 7 C A p Q 7 + 7 ( 2 A p Q$ + m 7 rj X p Q 9 + k ( rj X p Q w + 7 V X p 2 Q u ) 



a 3 = {irh o 2 r)XR + m o ClXR l + 7™oCCAi? 2 + A; ( 2 C A #3 + C 3 7 C A #4 + 7 C 2 C 2 A #5 

+k m Q ( r] X R 6 + ( 2 7 X p R 7 + 7 ( ( X p R 8 + 7 m r] X p R 9 + k(r]XpR 10 + ^1] p 2 XR n ) 



«3 = (7^o 2 ^Ai? + ^oC 2 7Ai?i+7^oCCAi? 2 + A;C 2 CAi?3 + C 3 7CA J R4 + 7C 2 C 2 A J R 5 

+k m a ( r] X R 6 + ( 2 7 X p R 7 + 7 ( ( X p R 8 + 7 m r] X p R 9 + k ( i] XpR 10 + ^r] p 2 XR U ) 
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& 4 = ^-(m 2 1 r ] 2 So + ( 2 l( 2 S 1 +^(( 3 S 2 + m ( 1 (r ] S 3 + ^m ( 2 r ] S 4 + k(( 2 r ] S 5 

+km (r] 2 S 6 + ClCvpSi + *f( 2 r}pS 8 + m -fr] 2 pS 9 + k(r} 2 pS 10 + 7r? 2 p 2 Sn) 

«4 = y^(7™oV£ + C 3 7CSi+7C 2 C^^ 

+km (i] 2 S & + */((r]pS 7 + C 2 'yi]pS 8 + ; ym i] 2 pS 9 + k(rj 2 pS 10 +7^ 2 p 2 £n) 

«5 = (C 2 7CAT + 7CC 2 ATi + m ('jr]\T2 + ; ym (r]\T 3 + k((r]\T 4 
+krfi r ] 2 \T 5 + ( 1 r ] \pT (i + ^(r ] \pT 7 + kr ] 2 \pT 8 ) 

«5 = (C 2 7CAf + 7CC 2 Af 1 + m o C7^Af 2 + 7m o Cr7Af 3 + A;CC^Af4 
+A;m r ? 2 Af 5 + C7^Apf 6 + 7Cr7Apf 7 + A;r ? 2 Apf 8 ) 

«6 = (C 2 7CA^ + 7CC 2 A^i + m o C7^Af/ 2 + 7m o Cr?Af/ 3 + A;CC^Af/4 + A;^o^ 2 Af/5 
+C7^Apf/ 6 + 7Cr?Apf/ 7 + A;r ? 2 Apf/ 8 ) 

«6 = (C 2 7CAt7 + 7CC 2 A£7i + m ( 1 ri\U 2 + im (ri\U 3 + k((r ] \U 4 + krfi ri 2 \U 5 
+C-fr]XpU 6 + ^Cr]XpU 7 + kr] 2 XpUs) 





= -1 + X 




*(2,1) 


= -3 + 4x + 3x 2 




*(2,2) 


= -1 + 5 a; 2 




^(2,3) 


= -l + 2x + x 2 




£(2,4) 


= 1 — 6x + 7 x 2 




£(2,5) 


= l-5x + 2x 2 




£(2,6) 


= -3 + 3x + 2x 2 




*(3,1) 


= —3 + x + 5 x 2 4 


-3 a: 3 


£(3,2) 


= (-l + 3x) (l + 


a: 2 ) 


£(3,3) 


= 2-9x + 6a; 2 + 5a; 3 


£(3,4) 


= -2 + 9 a; -9 a; 2 


+ 4x 
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t (4jl) = 3- I6a; + 21x 2 -6x 3 + 2x 4 

t (4j2) = \-x-5x 2 + 5 a; 3 + 4 x 4 

t (43) = 2-7x-3x 2 + llx 3 + 13x 4 

t (4)4) = i_5 ;r + 6x 2 -5x 3 + 7x 4 

* (5i i) = -l + 9x- 25x 2 + 29x 3 - 18x 4 + 14x 5 

t (5j2) = _2 + 8x-3x 2 -2x 3 -8x 4 + x 5 

t (5)3) = 3-3x-15x 2 -x 3 + 30x 4 + 14x 5 

t (5)4) = -l-3x + llx 2 + 8x 3 -12x 4 + 5x 5 

t (55) = _3 + g a; + 5 a ;2_ 10a ,3_ 14x 4 + a; 5 

t (6>1) = l + 4 ;r _38a; 2 + 47x 3 + 10x 4 - 23x 5 + 7x 6 

t (7>1) = l- 8x + 25x 2 - 44x 3 + 64x 4 - 74x 5 + 36x 6 + 12x 7 

t (7>2) = 3 - 21x + 60x 2 - 114x 3 + 181x 4 - 147x 5 + 4x 6 + 18x 7 

t (7>3) = _l + 5 a; -3 ;r 2 _ 2x 3 - 39x 4 + 69x 5 - 37x 6 + 40x 7 

t (7>4) = 3 - 22x + 61x 2 - 102x 3 + 157x 4 - 142x 5 + llx 6 + 18x 7 

t (7i5) = 2 + llx- 72x 2 + 34x 3 + 123x 4 - 30x 5 - 49x 6 + 13x 7 

t (8 ,i) = 2 - 35x + 183x 2 - 398x 3 + 474x 4 - 441x 5 + 387x 6 - 170x 7 + 54x 8 

t (8>2) = i_5 ;r + 9 ;r 2 _i0x 3 - 33x 4 + 195x 5 - 245x 6 + 44x 7 + 12x 8 

t (83) = 3 - 77 x + 528 x 2 - 1660 x 3 + 2967 x 4 - 3417 x 5 + 2702 x 6 - 1342 x 7 + 360 x 8 

t (8>4) = -l + 19x- 115x 2 + 286x 3 - 211x 4 - 341x 5 + 767x 6 - 548x 7 + 128x 8 

t (8>5) = -3 + 23x- 81x 2 + 191x 3 - 307x 4 + 261x 5 - 91x 6 + 21x 7 + 18x 8 

t (9>1) = _i + i2;r- 39x 2 + 30x 3 + 79x 4 - 172x 5 + 63x 6 + 50x 7 - 6x 8 + 48x 9 

t (9)2) = l- 16x + 81x 2 - 174x 3 + 118x 4 + 222x 5 - 529x 6 + 500x 7 - 303x 8 + 84x 9 

t (9)3) = -l + 21x- 135x 2 + 405x 3 - 599x 4 + 337x 5 + 127x 6 - 173x 7 - 128x 8 + 210x 9 

t (9>4) = i_21x + 97x 2 - 155 x 3 + 87x 4 - 29 x 5 + 3 x 6 + 63 x 7 - 108x 8 + 30x 9 

t (10)1) = 5 - 51 x + 219 x 2 - 537 x 3 + 883 x 4 - 1025 x 5 + 757 x 6 - 287 x 7 + 16 x 8 + 12 x 9 + 72 a: 

t (10>2) = -1 + 21 x - 152 x 2 + 572 a; 3 - 1335 x 4 + 2176 a; 5 - 2631 x 6 + 2340 x 7 
-1380 a; 8 + 315 a; 9 + 171a; 10 

i (10)3) = -3 + 28 a; - 150 x 2 + 486 x 3 - 8 10 a; 4 + 422 x 5 + 428 a; 6 - 670 x 7 
+605 a; 8 - 394 a; 9 + 186 a; 10 

t (10 4) = 1 - 16 x + 96 x 2 - 338 a; 3 + 864 x 4 - 1546 a; 5 + 1426 x 6 + 114 x 7 
-1377 a; 8 + 762 a; 9 + 270 x 10 

4(11,1) = -1 + 18 x - 121 x 2 + 482 a; 3 - 1294 a; 4 + 1958 a; 5 - 226 a; 6 - 4462 x 7 + 7239 a; 8 
-4456 x 9 + 483 x 10 + 252 x 11 

t (12)1) = 3 - 56 x + 383 x 2 - 1302 x 3 + 2460 x 4 - 2908 x 5 + 2826 x 6 - 2040 x 7 
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-2121 x 8 + 8532 x 9 - 9449 x 10 + 2766 x 11 + 1674 x 



P2 = {-l + 2x)(l+x) 

p 3 = _l + l0a;-17a; 2 + 12a; 3 

p 4 = (-1 + 3x) (-1 +5x + x 3 ) 

p 5 = l-7x + 21a; 2 -32a; 3 + 20a; 4 + 9x 5 

P = -12p 3 P5f(i,i) Pi = 24x 3 t (lil) t (10 ,i) i 3 2 = -24a;t ( io ) 2)t(i i i) 

P 3 = 4xt ( i lj i ) t 2 11) P 4 = 8p 3 p 5 t(2,3)*(i,i) ^5 = 4x 2 p 3 p 5 tJ lil) 

«i = (m 2 r ? 2 APo + C 2 C 2 APi + m CC^AP 2 + CC^ApP 3 + m r ? 2 ApP 4 + r ? 2 Ap 2 P 5 ) 

Qo = 6"s/3p 3 P4*(i,i) <5i = 6^3; (-1 + 3x) £( 8 ,i)£(i,i) 

Q 2 = 6^ x 2 (-1 + 3 x) p 2 2 t( 3l i)*(i,i) ^3 = -12 (-1 + 3 x) t (hl) i (9)1) 

—6 o 2 —6 o 9 2 

g 4 = X (-3 + 5x) t (3i2 )t(4,l)t(i i i) Q 5 = X {-3 + 5x) p 2 £(1,1) £(3,2) 

g 6 = -12>/3a; 2 (-1 + 3x) p 3 t(4,2) *(i,i) Q 7 = -2>/3a; (-l + 3x) * (9) 2)t(i,i) 

Q 8 = -6\/3x 3 (-l + 3i)p 2 2 i(2,6)4,i) = -4\/3p 3 P4 *(2,3) *(i,i) 

Q10 = 4v / 3x 3 (-1 + 3x) P3 *(i,i) *(5,3) Q11 = ~ 2 y/3x 2 p 3 p4, £( 1>1} 

Qo = 6V3p 2 P5^(i,i) Oi = 12 >/3a;p2*(7.i)*(i,i) 

(Q 2 = -12^^ (-l + 3x) p 2 t {5 , 2) tf hl) Q 3 = 12V3x 3 p 2 t {7 , 2) tl ltl) 

- —6 v^3 3 3 2 ^ —6 3 2 

V4 = X £(3,2) £>2 £(i,i) <^5 = X p 2 £(3,2) £(4,1) *(1,1) 

Oe = 24 V / 3x 2 p 2 p 5 £( 11 ) Q 7 = -2 v^x^ £(8,2) £(1,1) 

Qs = 2v / 3x 2 (-1 + 3x) p 2 £(6,i) £(1,1) Oq = -4 V^Ps £(2,3) £(1,1) 

O10 = -4\/3a; 2 (1 + x) P2P5 £(2,4) £(1,1) On = ~ 2 v / 3x 2 p 2 P5 £(1,1) 
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Ro = 6a/3p2P5*(i,i) 

R 2 = 12v / 3xp 2 i(7,i)^(i,i) 

—6 3 2 

*U = X p 2 £(3,2) *(4,1) *(i,i) 

# 6 = 24v / 3x 2 p 2 p 5 t(i,i ) 

i? 8 = -2 V3ip2%2)t(i 1 i) 

i?l = -4v / 3x 2 p 2 P5 (1 +X) t(2,4)*(i,i) 



i?l = -12 X 2 (-1 + 3 x) P 2 £(5,2) 

R 3 = 12VSx 3 p 2 t {7 , 2) £ 2 ljl} 

— 6 V^3 33 2 
^5 = ^ P2 ^(3,2) ^(1,1) 

i? 7 = 2v / 3a: 2 (-1 + 3 a;) P2 f( 6) i) 
R g = -4 v / 3p 2 P5i(2,3)i(i,i) 
i?n = -2"\/3a; 2 p2P5*(i,i) 



-Ro = 6v / 3p 3 P4^i,i) 

i? 2 = 6V3x (-1 + 3 a;) *(8,1) *(!,!) 

^4 = 6 ^ a; 3 (-3 + 5 x) p 2 t 2 {11) £( 3i2 ) 

Rq = -12VZX 2 (-1 + 3 a) p 3 *(4,2)*(i,i) 

i? 8 = -2 V / 3x(-l + 3 x) £( 9 ,2) *(i,i) 

i?10 = 4V3X 3 (-l + 3x) P3*(i,i)f(5,3) 



i?l =6V / 3X 2 (-1 + 3 a) P2 2 *(3,l) *(!,!) 



i? 3 = -12 (-1 + 3 x) i ( i,i) * (9 ,i) 

i? 5 



6 \/3 

X 3 (~3 + 5x) t(i,i) t( 3) 2) f(4,l) 



r/ 



i? 7 = -6 v^3 x 3 (-1 + 3 x) p 2 2 £ (2 , 6) tj ljl} 

#9 = -4^^3^4^(2,3)^(1,1) 

#n = -2 v / 3x 2 p 3 p 4 i(i,i) 



So = -6 v / 2p 3 o-i(3,4)i(i,i) 

S 3 = -3V2xat(s t s) £ 3 1;1) 

S e = 6 v / 2a;a-p 3 t(5,i)t(i i i) 

Sg = 4v / 2p 3 <7£( 2 ,3) £(3,4) £(l,l) 



5i = -18v^a: 3 ^(4,i) *(i,i) 
S 4 = -9\/2xap 2 2 t( 2 ,5)*(i,i) 
S 7 = 3 v / 2a;(3-t(8,4) *(i,i) 

5i = -2v^X^P3*(i,i)*(7,3) 



S 2 = -18v^^ 3 p2 2 t(i,i) 
S 5 = 12 v / 2a; 3 o-t(9 i3 ) t(i,i) 
S 8 = 3 v 7 ^ x a p 2 t (4i 4) t( ljl} 

5n = 2 X 2 O- p 3 £( 3)4 ) t 



S = -6V2ap 3 t (3A) tf 1A) 

S 3 = -3 v / 2xo-t(8, 3 )t 3 lil) 

S 6 = 6V2xap 3 t {5A) t 3 {11) 

Sg = 4V2&p 3 t( 2 ,3)t(3,4)t(l,l) 



Si = -18v / 2x 3 op 2 2 i(i,i) 

S 4 = -9 X (X p 2 2 £(2,5) *(i,i) 

S 7 = 3 V2xat (8A) tf 1A) 
Sio = -2 V2 xa p 3 t 2 hl) t( 7)3 ) 



S 2 = -18v^x 3 a* (4i i) 

S 5 = 12 V / 2x 3 0-t( 9i3 )t(i i i) 

S 8 = 3 V2 x a p 2 £ (4i 4) tj ljl} 
Sn = 2 \/2 x 2 a p 3 t( 3j4 ) £(1,1) 
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T = 6x 2 (-1 + 3x) t ( i,i) t(io,3) T 1 = 18 x 2 p 2 t {7A) tf ltl) 

T 2 = 6x (-1 + 3x) P3 ^(5,4) *(i,i) r 3 = 18a;p2P5<(i,i) 

T 4 = -12a; 2 t(i2,i) T 5 = -12 xp 3 p 5 t 2 ltl) t (2t2) 

T 6 = -2x (-1 + 3x) P3 *(i,l) *(7,5) r 7 = -6xp 2 P5^(3,3)i(i,i) 

T 8 = 4 xp 3 p 5 i( 4 ,3) 

T = 18 x 2 p 2 1( 7 ,4) Ti = 6 x 2 (-1 + 3 x) £(i , 3 ) 

f 2 = 18xp 2 P5*(l,l) ^3 = 6X (-1 + 3x) P3 *(5,4) 

T 4 = -12x 2 i(i 2i i) T 5 = -12 xp 3 p 5 t( 2j 2) t 2 !,!) 

f 6 = -6 Xp 2 p 5 t ( 3)3 ) f 7 = -2X (-1 + 3x) P3 t(7,5) *(1,1) 

T 8 = 4xp3p 5 t(i,i)t(4, 3 ) 

C/ = 6 V3ix 2 (-1 + 3x) t(9, 4 ) t 2 1;1) f^i = 6 v / 3zx 2 p 2 t( 8j5 ) t 2 1;1) 

C/ 2 = -6 V^i X (-1 +3 X) p 3 *(5,5) ^3 = -6\/3zxp2P5*(2,i)*(i,i) 

^4 = -12v / 3ix 2 t(i 0i 4)t 2 11) C/5 = 12 v3ia;p3P5*(i,i) 

C/ 6 = 2 v^ix 2 (-1 + 3x) p 3 t( 5i 4)t( li i) C/ 7 = 6 V3ix 2 p 2 p 5 tf hl) 

U s = -4v / 3ix 2 p 3 p 5 t( 2i2 )t 2 lil) 

Uq = -6 v / 3ix 2 p 2 t( 8 ,5)i(i i i) f/i = -6 v^ix 2 (-1 + 3x) t( 9j4 ) t 2 l5l) 

U 2 = 6V3ixp 2 p 5 t {2:1) tf 11) U 3 = 6\/3ix (-l + 3x) Pa *(i,i) *(5,b) 

t/ 4 = 12 v / 3«x 2 t(io, 4 ) i(i,i) C/ 5 = -12 \/3ia;p3P5*(i,i) 

C/e = -6v / 3zx 2 p 2 p 5 t( 11) t7 7 = -2 v^ix 2 (-1 + 3x) p 3 t( 5i 4) ^i jl} 

C/ 8 = 4v / 3ix 2 p 3 p 5 t( 2i2 )t 2 1)1) 
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